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INTRODUCCION  I 


ANTECEDENTES  HISTORICOS 


L interés  por  el  procesado  digital  de  imágenes  se  divide 
en  dos  áreas  principales  de  aplicación:  mejora  de  imagr- 
•¿i  nes  (para  su  mejor  interpretación  por  el  hombre)  y pru- 

■ll^B  cesado  de  imágenes  (para  su  interpretación  automalica 
por  medio  de  una  máquina).  Las  primeras  técnicas  de  tra- 
tamiento de  imágenes  se  remontan  al  principio  de  los  años 
veinte,  cuando  se  trataba  de  mejorar  las  fotograbas  de  pe- 
riódicos enviadas  por  cable  submarino  entre  Londres  y 
Nueva  York,  El  uso  de  imágenes  en  forma  digital  se  re- 
monta a 1921.  con  la  introducción  del  primer  sistema  de  transmisión  de 
fotografías  por  cable  submarino  desde  Londres  a Nueva  York.  La  miro  uc- 
eión  del  sistema  Barlanc  de  transmisión  fotográfica  por  cable  redujo  el 
tiempo  empicado  en  transportar  una  fotografía  a través  del  Atlántico  de 
una  semana  a menos  de  tres  horas.  Los  primeros  sistemas  Bariane  eran  ca- 
paces de  codificar  imágenes  en  cinco  niveles  distintos  de  intensidad  y esta 
capacidad  fue  aumentada  en  1929  a quince  niveles.  Durante  este  periodo 
también  ve  mejoró  notablemente  d proceso  de  reproducción,  que  paso  de 
ser  fo  impresora  tclcgr^ficD  a un  sisl¿ma  con  placa  fotográfica  tmpresio 
nada  con  rayos  de  luí  modulados  por  la  imagen  codificada. 


Aunque  continuaron  las  mejoras  en  los  métodos  de  procesado  durante 
los  siguientes  treinta  y cinco  años,  la  gran  potencia  de  cálculo  de  los  com- 
putadores. derivada  de  la  progresiva  evolución  de  las  tecnologías  VLM  y 
el  programa  espacial  fueron  Sos  que  dieron  gran  impulso  a los  conceptos 
de  procesado  de  imagen.  El  tratamiento  digital  de  imágenes  entendido 
como  la  aplicación  de  procedimientos  matemáticos  tendentes  a corregir 
las  fotografías  y mejorar  su  interp-retábilídad  no  se  inició  hasta  1963.  cuan- 
do en  los  laboratorios  de  investigación  de  IBM  se  empezó  a procesar  me- 
díame ordenador  las  imágenes  obtenidas  por  el  satélite  Tiros’  A partir  dé 
ese  momento,  el  tratamiento  digital  de  imágenes  se  ha  convertido  en  un 


5 


área  de  las  ciencias  aplicadas  que  viene  experimentando  un  crecimiento 
acelerado  duran  le  los  últimos  veinte  anos.  El  uso  de  técnicas  de  mejora 
de  imagen  por  ordenador  comenzó  en  el  Jei  Propulsión  Laboratory  [Pasa* 
dena.  California)  en  1964,  cuando  fotografías  de  la  Luna  transmitidas  por 
el  /fangcr  7 eran  procesadas  por  un  computador  para  corregir  varios  tipos 
de  distorsión  de  la  imagen  inherente  a la  cámara  de  TV  de  a bordo.  Estas 
técnicas  sirvieron  corrió  base  para  posteriores  métodos  usados  en  la  me- 
jora y restauración  de  imágenes  para  programas  tan  familiares  como  las 
misiones  Suryeyor  a la  Luna,  la  serie  de  misiones  Mariiicr  a Marte  y los  vue- 
los tripulados  Apollo  a la  Luna.  Este  desarrollo,  inicialmente,  se  efectuó 
con  objeto  de  tratar  las  imágenes  suministrada»  por  las  diversas  misiones 
espaciales,  en  especial  las  de  los  satélites  dedicados  a la  teledetección  de 
recursos  naturales  c impulsada  por  la  utilización  de  computadores  cada  vez 
más  potentes,  tuvo  menor  incidencia  en  otras  áreas,  si  bien  en  algunos  ca- 
sos, como  el  de  ja  microscopía  electrónica,  estas  técnicas  se  han  emplea- 
do desde  mediados  de  k década  de  los  sesenta. 

Desde  1964  hasta  ahora,  el  campo  del  procesado  de  imágenes  ha  expe- 
rimentado un  vigoroso  crecimiento.  Además  del  programa  espacial,  las 
técnicas  más  usadas  hoy  se  dedican  a una  gran  variedad  de  campos  como: 
medicina  (imágenes  biomedícas  y rayos  X),  geografía  y meteorología  (ana- 
tisis  de  imágenes  aereas  y de  satélite),  arqueología,  astronomía,  biología, 
medicina  nuclear,  defensa  o industria.  Todas  ellas  no  sólo  en  cuanto  a la 
mejora  de  imágenes  para  interpretación  humana,  sino  también  en  lo  refe- 
rente a su  interpretación  por  una  máquina  en  campos  como:  reconoci- 
miento de  Caracteres,  robots  industriales  para  montaje  e inspección,  reco- 
nocimiento de  objetivos  militares,  procesado  de  huellas  dactilares,  análi- 
sis de  muestras  de  sangre  y radiográficas,  predicción  del  tiempo  atmosfé- 
rico y de  las  cosechas,  acústica,  comunicaciones,  geofísica,  procesamien- 
to de  voz,  microscopía,  teledetección,  etc. 

Por  otra  parte,  se  ha  producido  un  gran  acercamiento  entre  las  técni- 
cas de  procesado  digital  de  imágenes  y la  de  gráficos  por  ordenador  /rom- 
puier^phicsl  hasta  d punto  de  que  ambas  se  solapan  en  el  marco  de  cier- 
tas aplicaciones. 

El  tratamiento  digital  de  imágenes  abarca  una  gran  variedad  de  técni- 
cas y herramientas  mate m aticas  que  se  han  desarrollado,  si  bien  cabe  se- 
ñalar que  en  los  últimos  años  se  han  comenzado  a desarrollar  bases  cien- 
tíficas  en.  lo  que  antes  era  un  simple  conjunto  de  técnicas. 

El  tratamiento  digital  de  imágenes  abarca  una  gran  variedad  de  técni- 
cas y herramientas  matemáticas  que  pueden  ser  clasificadas  en  los  siguien- 
tes puntos: 

I Digital ización  y cuanlifícación,  que  consiste  en  efectuar  el  paso  de 
imágenes  continuas  a discretas. 


% Operaciones  de  preprocesado,  consistentes  en  la  realización  de  de- 
terminadas operaciones  tendentes  a normalizar  los  datos  de  entrada. 

1 Mejora  o realce  de  imágenes,  consistente  en  efectuar  algunas  ope- 
raciones tendentes  a mejorar  la  apariencia  de  los  datos  de  entrada. 

4 Restauración  de  imágenes,  consistente  en  la  recuperación  de  la 
imagen  original,  en  el  supuesto  de  que  los  datos  de  entrada  se  encuentran 
afectados  de  algún  tipo  de  degradación,  como,  por  ejemplo,  desenfoque, 
turbulencia  atmosférica*  movimiento,  etc, 

5 Codificación  de  imágenes,  consistente  en  la  compresión  de  la  in- 
formación con  objeto  de  reducir  los  costes  de  transmisión  de  las  imáge- 
nes, asi  como  conseguir  un  almacenamiento  más  eficiente,  como  sucede, 
por  ejemplo,  en  el  caso  de  las  imágenes  médicas,  en  ingeniería,  etc. 

Ó Análisis  de  imágenes,  consistente  en  la  obtención  de  descripciones 
numéricas  o simbólicas  de  la  imagen  de  partida.  Las  descripciones  de  una 
imagen  se  refieren  generalmente  a las  panes  constituyentes  de  dicha  ima- 
gen,  Es  por  ello  que.  el  primer  paso  en  todo  proceso  de  análisis  de  imáge- 
nes So  constituye  la  denominada  segmentación  en  regiones  con  significa- 
do propio. 

7 Reconstrucción  de  imágenes,  consistente  en  la  obtención  de  infor- 
mación acerca  de  un  objeto  a partir  de  sus  proyecciones,  como,  por  eje  ni- 
pío.  la  realizada  en  tomografía,  en  donde  se  obtienen  secciones  transver 
sales  del  cuerpo  humano,  o los  métodos  basados  en  la  resonancia  magné- 
tica nuclear. 

Estas  técnicas  adoptan  en  la  práctica  formas  y modos  de  operación  muy 
diversos  en  función  del  tipo  de  imágenes  a que  se  aplican  y del  tipo  de  in- 
formación que  de  ellas  se  desea  extraer. 


REPRESENTACION 
DIGITAL  DE  IMAGENES 

Antes  de  comenzar  con  fa  descripción  de  algunas  de  las  técnicas  del 
procesado  digil  al  de  imágenes,  es  preciso  definir  que  se  entiende  por  una 
imagen  digital,  a fin  de  poder  apreciar  las  posibilidades  de  estas  técnicas, 
asi  como  comprender  la  problemática  que  entraña  su  utilización. 

En  el  ámbito  del  iraiamiertto  de  imágenes,  el  término  im&g*n  se  refiere 
a una  función  b ¡dimensional  de  la  intensidad  luminosa  f(x,y).  donde  x e v 
denotan  coordenadas  espaciales  y el  valor  de  f en  cada  punto  (x,y)  es  pro- 
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poiuonal  al  bríUo,  intensidad  □ nivel  de  gris  de  la  imagen  en  esc  punió,  ’ 
A vetes  es  útil  visualizar  una  función  imagen  en  perspectiva  iri dimensio- 
nal, siendo  la  intensidad  el  tercer  eje. 

Una  bnagen  digital  f (x,y)  se  obtiene  a partir  de  la  radiancia  de  un  ob- 
jeto que  ha  sido  diserer izada  en  ambas  coordenadas  espaciales  y en  míen-  ¡ 
stdad  Podemos  considerar  una  imagen  digital  como  una  matriz  cuyos  ín- 
dices de  fila  y columna  identifican  un  punto  de  la  imagen  y e]  correspon- 
diente elemento  de  la  nutria  representa  el  nivel  de  gris  de  ese  punto.  Los 
elementos  de  dicha  matriz  digital  se  llaman  elementos  de  imagen  o pixeh 
<o  pcls),  como  abreviatura  de  *pieturc  clements». 

Aunque  el  tamaño  de  una  imagen  digital  varia  con  la  aplicación,  hay 
grandes  ventajas  en  elegir  matrices  cuadradas  cuyos  tamaños  y número  de 
niveles  de  gris  sean  potencias  enteras  de  2.  Por  ejemplo,  un  tamaño  típico 
comparable  en  calidad  a una  mugen  de  TV  monocromática  es  el  de  una 
matriz  de  51 2 x 512  con  128  niveles  de  gris. 

Una  de  las  características  que  define  el  procesado  de  imágenes  es  el  ele- 
vado volumen  de  información  que  conlleva  este  tipo  de  procesos.  Este  he- 
cho hace  que  el  coste  compuiadonal  de  las  operaciones  que  se  efectúen 
con  imágenes  sea  muy  elevado,  a menos  que  éstas  sean  muy  sencillas. 


ELEMENTOS  DE  UN  SISTEMA  DE  PROCESADO 
DIGITAL  DE  IMAGENES 

Los  elementos  de  un  sistema  de  procesado  de  imagen  básico  v de  pro- 
pósito general  son  los  de  la  figura  1. 


A I macenaniie  nto 


Imagen 

- 

Digitalizador 

Ordenador 

Digital 

- 

Consola  de  , 
Operación 

Visual  izador 


t-ttf.  I fibnmum  de  un  tHxttmo  de  procesamiento  digital  de  imágenes. 
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Ij.  operación  de  un  sistema  de  este  tipo  puede  dividirse  en  tres  catego- 
rías principales:  digilalizaeiún,  procesado  y visual  ilación. 


a|  Digitalizadartts 

Convierten  una  imagen  en  su  representación  numérica,  de  forma  que 
pueda  ser  introducida  en  un  ordenador  digital,  Los  dispositivos  digitaltza- 
dore*  más  comunes  son: 

Mlc  rodé  nal  lómeteos.  En  ellos  la  transparencia  o fotografía  se 
monta  en  una  plataforma  sobre  la  que  se  lanza  un  fino  ha/  de  luz.  La  plá* 
i .llorína  se  hace  girar  o mover  y se  van  obteniendo  los  valores  puntuales. 
E¡  rayo  reflejado  (o  el  que  la  atraviesa,  en  el  caso  de  transparencias)  es  en- 
ligado hacia  un  fotodelector  y el  nivel  de  gris  es  codificado  en  base  a la 
intensidad  del  rayo  en  cada  punto.  Una  imagen  digital  se  obtiene  permi- 
tiendo sólo  valores  discretos  de  intensidad  y posición  a la  salida  Son  dis- 
positivos  lentos  pero  de  ütxa  gran  precisión. 

S canoera  «Flylng  Spot»,  También  operan  sobre  el  principio  de 
enfocar  un  haz  de  luz  reflejado  o transmitido  sobre  un  foto-detector.  En 
este  caso,  sin  embargo,  la  imagen  permanece  estática  y la  fuente  luminosa 
es  un  tubo  de  rayos  catódicos  (CRT)  en  el  cual  un  haz  de  electrones  incide 
sobre  una  superficie  fluorescente.  El  rayo  produce  un  punto  de  luz  que  se 
mueve  y analiza  la  imagen  sobre  la  superficie  del  tubo.  Ya  que  el  rayo  se 
mueve  electrónicamente,  permite  grandes  velocidades  de  lectura.  Tam- 
bién permite  gran  flexibilidad  a la  hora  de  leer  una  imagen,  ya  que  la  po- 
sición de Í rayo  puede  ser  rápida  y fácilmente  establecida  extemamente. 

Disee torea  de  Imagen,  Su  operación  se  basa  en  el  principio  de 
emisión  electrónica. 

— Cámara»  de  TV.  Muchos  dígita  liza  dores  de  cámaras  de  TV  em- 
pican un  tubo  vidicón,  cuya  operación  se  basa  en  el  principio  de  la  foto- 
conduetlviidád.  Una  imagen  enfocada  en  la  superficie  del  lubo  produce  un 
patrón  de  variación  de  conductividad  que  corresponde  a la  distribución 
de  intensidad  de  la  imagen  óptica.  Éstos  dispositivos  son  menos  precisos 
que  los  anteriores,  pero  generalmente  más  baratos. 


b)  Procesa  dores  de  imagen 

Los  dispositivos  usados  para  el  procesado  de  imagen  varían  desde  mi- 
croprocesadores para  aplicaciones  de  propósito  general  hasta  enormes  sis- 
temas de  computadores  capaces  de  realizar  gran  variedad  de  fuñe  iones  en 
imágenes  matrieiales  de  alta  resolución  (coni  pul  adores  de  propósito  espe- 
cial). 


c)  Visittllizatíurvi 

Son  los  dispositivos  que  se  encardan  di  convenir  las  matrices  numéri- 
cas almacenadas  en  el  ordenador  en  una  representación  adecuada  para  la 
interpretación  humana-  Hay  tres  tipos  principales  de  dispositivos  vi  su  aliga- 
dores: 

— C.R.T,  Los  tubos  de  rayos  catódicos  (CR.T.)  son  sistemas  en  los 
que  la  posición  horizontal  y vertical  de  cada  elemento  de  la  matriz  imagen 
es  convertida  en  voltajes  que  se  usan  para  dirigir  un  haz  de  electrones,  de 
forma  que  se  produce  Ea  imagen  de  salida  b [dimensional,  En  cada  punto 
de  Ea  pantalla,  la  intensidad  del  haz  es  modulada  por  un  voltaje  proporcio- 
nal at  valor  del  punto  correspondiente  en  la  matriz  (pixcl). 

— Monitores  de  TV,  Estos  dispositivos  convierten  la  imagen  alma- 
cenada en  el  ordenador  en  una  señal  de  vídeo  que  puede  ser  visualizada. 
Su  ventaja  es  que  pueden  ofrecer  una  gama  de  tonalidades  que  se  aproxi- 
ma en  calidad  a la  de  una  placa  fotográfica,  produciendo  una  salida  que 
puede  ser  fácilmente  asimilada  por  el  sistema  visual.  Su  inconveniente  es 
que  el  monitor  debe  tener  un  ciclo  de  refresco  de  la  imagen  con  una  fre- 
cuencia de  unos  30  « brames»  (marcos  o imágenes)  por  segundo  para  evi- 
tar el  parpadeo»  mientras  que  la  mayoría  de  los  ordenadores  convencio- 
nales no  son  capaces  de  transferir  datos  a esta  velocidad-  Para  solucionar 
esto»  se  dota  al  sistema  de  TV  de  un  medio  de  almacenamiento  o «buHer» 
para  transferir  datos  al  monitor  a la  frecuencia  de  video. 

— Impresoras*  Los  dispositivos  de  impresión  de  imagen  son  útiles 
para  trabajos  de  procesado  de  imágenes  a baja  resolución.  Un  medio  de  ge- 
nerar los  diferentes  niveles  de  gris  de  la  imagen  directamente  en  el  papel 
es  aprovechando  la  capacidad  de  *over-slrike»  (sobre-impresión)  de  las  im- 
presoras convencionales.  El  nivel  de  gris  de  cada  pixel  puede  ser  contro- 
lado por  el  número  y Ea  intensidad  de  los  caracteres  sobre  impresionados 
en  él.  Seleccionando  adecuadamente  el  conjunto  de  caracteres,  puede  con- 
seguirse una  distribución  de  grises  aceptable  con  un  programa  sencillo  y 
relativamente  pocos  caracteres.  Este  sistema  de  sobreimpresión  de  carae 
!ere$  es  el  que  se  ha  empicado  para  visualizar  Eas  imágenes  definidas  en 
ios  programas  del  apéndice  2,  proporcionando  un  grado  de  resolución  su- 
ficiente para  el  propósito  de  entrenamiento  al  que  se  hallan  dirigidos. 


to 


OPERADORES  SOBRE  IMAGENES 


MUESTREO 

Y CUANTIFIC  ACION  DE  IMAGENES 

ARA  poder  procesar  una  imagen  f(xry),  esta  función  debe 
pasar  por  un  proceso  previo  de  digitalización  tanto  espa- 
cial como  en  amplitud.  La  digital  ízación  de  las  coordena- 
das espaciales  es  lo  que  se  llama  maestreo  de  la  imagen 
y la  digital ización  en  amplitud  (cuantiheación  en  niveles 
de  grises) 

Una  imagen  continua  la  podemos  aproximar  por  una 
matriz  de  dimensiones  N x Ñ.  de  forma  que: 

pF(G,Q)  f<0(!)  * ftfLN-l)  i 

Í(\JQ)  F(  1 .1)  - f(LN-1) 

[fíN-1,0)  F(N-U)  .«  fíN-LH-1)  i 

En  esta  matriz,  cada  elemento  f(x»y)  es  un  valor  discreto.  Una  matriz 
de  este  tipo  se  denomina  imagen  digital,  y cada  uno  de  los  elementos  re- 
cibe el  nombre  de  pixeh 

Para  realizar  el  proceso  de  digitalizadón,  hay  que  fijar  previamente  el 
valor  de  N,  así  como  el  número  de  niveles  de  grises  que  se  van  a permitir 
para  cada  píxeL  Normalmente»  se  eligen  para  estos  valores  potencias  en- 
teras de  dos.  de  manera  que:  N-2"  y G-2",  siendo  G ct  número  de  niveles 
de  gris.  Según  esto,  el  número  de  bits  necesario  para  almacenar  una  ima- 
gen digitalizada  será  de  N x N pixels  de  dimensión  y m bits/ptxel  será: 
b - N x N x m 

Así.  por  ejemplo,  para  almacenar  una  imagen  de  128  x 128  pixels  y de 
64  niveles  de  gris  se  necesitarán  98.3Ü4  bits. 

Un  factor  importante  a tener  en  cuenta  es  la  determinación  de  cuántos 
niveles  de  gris  y qué  valor  de  N se  necesitan  para  obtener  una  buena  aproxi 
rriación  de  la  imagen.  La  resolución  depende  mucho  de  los  valores  que  se 
den  a m y a Nr  de  forma  que  cuanto  más  altos  sean  éstos,  mayor  sera  la 
resolución  de  la  imagen.  Sin  embargo,  los  requerimientos  de  memoria  y 


f(xry)- 
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proceso  numen  tan  considerablemente  aJ  aumentar  estos  parámetros.  Para 
obtener  una  imagen  de  calidad  análoga  a la  de  TV  se  necesitan  512x512 
pifcels  y I2S  niveles  de  grises.  Como  norma  general,  un  sistema  que  pro- 
cese imágenes  debe  ser  capaz  de  mostrar  imágenes  de  256  x 256  pixels  y 
64  niveles  de  grises,  como  mínimo. 

El  efecto  que  produce  fa  reducción  de  N es  una  imagen  digital,  es  la 
aparición  de  un  efecto  damero  en  Ea  misma.  Al  disminuir  m,  el  efecto  que 
se  produce  es  la  aparición  de  falsos  contornos.  Estos  son  los  efectos  que 
produce  la  variación  de  N y m por  separado;  en  1965  Huang  estudió  las 
relaciones  entre  los  dos  parámetros,  hadándolos  variar  aun  tiempo.  Para 
este  experimento  utilizó  tres  imágenes  con  distintos  niveles  de  detalle:  a 
las  conclusiones  a que  llegó  fueron: 

L Al  aumentar  N y m aumenta  la  calidad  de  imagen. 

2.  En  algunos  casos,  para  un  N fijo  se  mejora  la  imagen  al  disminuir 
m,  esto  es  debido  a que  se  aumenta  el  contraste. 

3.  En  imágenes  con  una  gran  cantidad  de  detalle,  se  necesitan  sólo 
unos  pocos  niveles  de  grises  para  conseguir  una  buena  aproximación. 


- MUESTREO  Y CUANTIFICACION  NO  UNIFORMES 

Fijado  el  valor  de  N,  en  muchos  casos  es  posible  mejorar  el  aspecto  de 
Ja  imagen  usando  un  método  adaptativo.  en  el  que  el  proceso  de  muestreo 
depende  de  las  características  de  la  imagen.  En  general,  se  necesita  un  ma- 
yor número  de  muestras  en  aquellas  regiones  en  que  las  transiciones  en- 
tre distintos  niveles  de  grises  son  más  pronunciadas,  haciéndose,  por  tan- 
to, un  muestreo  más  preciso,  pudiéndose  tomar  un  número  menor  de 
muestras  en  las  regiones  de  la  imagen  donde  las  transiciones  son  más 
suaves. 

Vamos  a considerar,  por  ejemplo,  la  imagen  de  un  rastro  sobre  un  fon- 
do uniforme.  El  fondo  queda  representado  adecuadamente  si  tomamos  un 
número  pequeño  de  muestras,  mientras  que  no  ocurre  lo  mismo  con  el  ros- 
tro, por  lo  que  hay  que  aumentar  el  número  de  muestras  tomadas  de  este 
ultimo,  mejorándose  asi  la  Imagen. 

Se  ve  claramente  que  hay  imágenes  en  las  que  distribuir  el  número  de 
muestras  de  forma  no  uniforme  resulta  ventajoso.  Una  buena  técnica  con- 
siste en  tomar  un  número  mayor  de  mu  es  Iras  en  Jas  fronteras  de  transi- 
ción entre  niveles  de  grises. 

El  muestreo  no  uniforme  se  usa  para  resaltar  los  bordes:  es  por  esto 
que  no  debe  usarse  en  imágenes  con  regiones  uniformes  relativamente  pe- 
queñas, por  ejemplo,  la  imagen  de  una  masa  de  gente  apiñada.  También 
se  usa  esta  técnica  cuando  se  quiere  obtener  una  imagen  más  detallada. 


Se  usa  uhh  técnica  parecida  a la  expuesta  anteriormente  en  el  mues- 
litti  no  uniforme,  pero  jugando  aquí  con  la  distribución  de  los  niveles  de 
grises. 

Hay  varias  formas  de  realizar  una  cuaniificaeión  no  uniforme  sobre  una 
imagen.  Una  de  ellas  consiste  en  utilizar  pocos  niveles  en  los  bordes,  de 
manera  que  se  consiga  resaltarlos,  y tomar  más  niveles  en  los  cambios  sua- 
ves. consiguiendo  la  reducción  y en  algunos  casos  ta  desaparición  de  fal- 
los contornos. 

Otra  alternativa  consiste  en  realizar  una  mayor  cuantificadón  en  bs  ni- 
veles de  grises  que  sean  más  frecuentes,  haciendo  una  cuanliRcación  me 
ñor  en  aquellos  niveles  que  sean  menos  frecuentes, 

Al  igual  que  el  muestreo  no  uniforme,  la  cuanlificaeión  no  uniforme  se 
usa  para  resaltar  bordes  y obtener  imágenes  con  mayor  detalle  de  represen- 
tación. 


LA  TRANSFORMADA  DE  FOURIER 

Uno  de  los  problemas  que  plantea  el  estudio  de  los  fenómenos  físicos 
es  el  de  su  representación  Dicha  representación  requiere  la  incorporación 
de  conceptos  que  son  idealizaciones  introducidas  por  el  observador  como 
mude  i ización  de  la  realidad.  Las  representaciones  más  habituales  de  una 
mal  son  la  temporal  (o  espacial)  y la  de  Frecuencias  temporales  (o  espa- 
i mies).  Esta  última  transformada  de  Fourier  de  la  primera,  permite  acce- 
der a las  componentes  espectrales  de  la  señal  (representación  armónica}. 
Ambos  tipos  de  representación  son  completas  y totalmente  equivalentes. 

Como  es  bien  conocido,  una  función  periódica  puede  ser  descompues- 
ta mediante  suma  de  infinitos  senos  y cosenos  con  diferentes  amplitudes, 
h rjLuencias  y fases.  En  el  caso  de  funciones  no  periódicas,  la  serie  infinita 
se  convierte  en  una  integral  de  la  forma: 

F(u)  - £ f(x)  exp(-j2rcux)  dx  (3.1} 

en  donde  exp(-]2xux)  - cm{2mix)  -jscn(2xuxh  siendo  j - I la  unidad 
imaginaria. 

La  Ec.  (3.1)  proporciona  la  transformada  de  Fourier  unidimensional 
F(u)  de  La  señal  f(x).  Si  ta  función  F(x)  es  real,  su  correspondiente  transfor- 
mada de  Fourier  es  compleja,  es  decir: 

(3-2) 


Ru)-R(u>+jI(u) 


3 


Fig.  2. 


Mádulo  de  tu  transformaéti  de  Fourier  ib)  de  imágenes  texturadxts  (a). 
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siendo  R(u)  e I(u)  tas  partes  rejal  e imaginaria  de  F(u).  También  puede  re- 
presentarse F(uJ  en  la  forma  módulo-arrumen  tal,  es  decir 

F(u)  - M(u>  exp(j0< u))  (hl) 

siendo  M(u)  = VR3  + I? 

y - tr>ci(ii)/R(u» 

Al  módulo  de  la  transformada  die  Fourier  F(u)  se  le  denomina  espectro 
de  Fourier  y a Oíu)  ángulo  de  fase.  En  la  figura  l se  presentan  algunos 
ejemplos  de  señales  con  variación  sinusoidal  y rectangular  (izquierda)  y 
su  correspondiente  transformación  de  Fourier  (derecha).  En  la  figura  2 se 
representan  los  módulos  de  la  transformación  de  Fourier  de  algunos  ejem- 
plos de  imágenes  con  información  de  textura. 


Transformada  de  Fourier  discreta 

Sea  ffx)  una  función  continua  que  se  ha  discretizado  en  la  secuencia 
iFta),  ffjtfl,  + 6),  f(x0  + 26) f(xo  + fN-l]$)|i,  tomando  N muestras,  separa- 

das por  6 unidades. 

Definimos  ahora  f(x)  - f(xfl+-  xfi),  y la  variable  x tomará  valores  discre- 
tos es  decir,  consideramos  la  secuencia  [F(Q),  f{l) f(N  1)1 

para  representar  los  valores  de  N muestras  equidistan  ciadas,  tomadas  de 
una  función  continua  f(x). 

Con  esta  notación,  la  transformada  de  Fourier  discreta  se  define  como; 

N-t 

F(u)  - 1/N  X f(?í)  cxp[-j2JCux/N], 

i-e 

para  u - 0,  1,2,  „„  N-  i , y: 


tf-3 

f(x)  » X F(u)  expfjZtaix/N], 

para  x ■ 0f  l.  2,  ,ut  N-  L 

Los  valores  u ■ 0,  1,  N-I  de  la  transformada  continua  se  correspon- 
den con  las  muestras  de  la  discreta  para  0,  $„  26,  (N-I)fi.  Por  otro  lado, 
la  relación  que  liga  S(x)  con  6(u)  viene  dada  por  la  siguiente  expresión; 

6(u)  - 1/N  5(x) 
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Extensión  caso  b [dimensional 


Para  dos  variables  la  transformada  de  Fourier  discreta  c$; 

N-  L K-1 

F(u,y)  - 1/MN  X £ £(x,y)  expl>j2it(ux/M  + vy/N)] 

■ *G  y-IJ 

pora  u - 0,  I. ...,  M-l  y v ■ 0,  1,  N~l 

M I H-l 

y r(x,y)  = X Z Hvt.v)  exp[)2x(ux/M  + vy/N>] 

u*tk  v-fl 

p.tra  x - 0,  1.  M-l  e y - 0,  1, ....  N-l 

Igual  que  en  el  caso  unidimensional,  f(x,y)  representa  muestras  de  la 
lu lición  f(x*  +■  xfix,  y9  + yby)  para  x - 0,  I,  ....  M-l  e y - 0,  1 N-l  Aná- 

logamente también  las  relaciones  entre  los  incrementos  en  los  dominios 
del  espacio  y la  frecuencia  son: 

6u«  1/M5* 
fiv  ■ I /N6y 

t liando  se  muestre*  una  señal,  se  obtiene  una  matriz  cuadrada;  por  tan- 
to, M - N,  y entonces: 

n i ht-i 

F(u,v)  - 1/N  X Xr<*y)  «ep[-j2*(i«  + vj)/N] 

y-fl 


para  u,v  » 0,  1,  N-l*  y. 


f(x.y)-  1/N  2 


Xf(u.v)  cxp(j2ir(ux  + vy}/N] 


pnra  x,y  - 0,  i N-l 

Hay  que  observar  en  estas  operaciones  que  se  ha  incluido  el  término 
1/N,  puesto  que  es  una  constante  multiplicativa  que  no  influye  en  el  par 
Hu.v)  y F(x,y)  transformadas. 
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Propiedad  de  separabílidad 

La  i rans  formada  discreta  de  Founer  para,  el  caso  bidimensional.  se  pur- 
de  expresar  como: 

N-|  N-l 

F{u,v)  “ l/K  2 exp  <-j2*ux/N)  X fíít.y)  cxp(-j2nvy/N) 

para  u,v  = 0,  1, N-1*  y: 

i*-t  »-i 

f{x,y)  • l/N  X CSP  (j2nux/N)  X f(u,v)  cxp(¿2rcvy/N) 
para  x.y  * 0,  J„  ,TÍ,  N-l 

La  ventaja  de  esta  propiedad,  es  que  F(u,v)  o ftx.y)  se  pueden  obtener 
en  dos  pasos  por  aplicación  sucesiva  de  la  transformada  de  Fourier  unidi- 
mensional o su  inversa.  Esto  se  ve  más  claro  si  hacemos  la  siguiente  sepa- 
ración de  términos: 

M-t 

F(u,v)  - 1/N  X F(x,v)  expH2rtux/N] 
n-fl 

en  donde. 

R*,v)  - N[l/N  £ f(x.y}  rxp(-]2jtvy/N)] 

En  esta  expresión,  para  cada  valor  de  x,  lo  que  nos  encontramos  entre 
corchetes  es  una  transformada  de  Fourier  unidimensional,  con  v - 0F  1 1 
N-l.  Por  tanto,  la  función  bidimensional  F(x,v)  se  obtiene  haciendo  una 
transformada  unidimensional  por  filas  sobre  f(x Fy)  y multiplicando  el  re- 
sultado por  N;  por  último,  para  obtener  Ffu.vj.  se  realiza  otra  transforma- 
da Unidimensional  por  columnas  sobre  F(x,v), 

Nótese  que  se  oble nd lian  los  mismos  resultados  haciendo  primero  le 
transformada  por  columnas  sobre  ffx.y)  y a continuación  haciéndola  por 
filas  sobre  el  resultado. 

Representación  de  la  transformada  de  Fourier 

Los  espectros  de  muchas  imágenes  decrecen  al  aumentar  la  frecuen- 
cia, por  lo  que  sus  términos  para  frecuencias  altas  se  hacen  impercepti- 
bles, Para  compensar  esto,  y lograr  una  mejor  visualiza*:  ¡ón,  se  utiliza  una 
representación  logarítmica  del  tipo: 

D(u,v)  - logfl  + K ab$(F(u,v))] 
en  lugar  de  Fíu.vJ.  La  función  D(u,v)  no  es  nunca  negativa. 


Con  esto  se  pretende  que  en  el  espectro  de  la  función  que  representa 
una  imagen  todos  los  picos  tengan  una  amplitud  similar,  de  forma  que  se 
cutí  siga  el  mismo  brillo  e intensidad  en  su  representaciún.Vamos  a demos- 
liarlo  viendo  que,  en  efecto,  la  relación  entre  d máximo  y el  mínimo  va- 
lor que  puede  tomar  la  función  es  menor  en  D(u,v)  que  en  F(u.v).  Sea: 

R - log(l  + F^J/logO  + FflUJ 

donde  F„,„  es  el  valor  máximo  alcanzado  por  Fíu.v),  y Fm „ es  el  mínimo 
de  tos  valores  que  alcanza. 

La  relación  F^/F^  es  siempre  mayor  que  R.  Por  ejemplo,  sea 
í ,,„  - 1,0,  y F^  - 0.2.  tenemos  entonces  F /F_._  - 5.0,  y R - 3.8;  si  ade- 

mas la  función  F(u,v)  la  multiplicamos  por  una  constante  k para  aumen- 
tarla, la  relación  F^^/F^  sigue  valiendo  5-0,  mientras  que  R,  si  k - 20. 
toma  el  valor  1.9.  que  es  mucho  menor, 


Transformada  rápida  de  Fourier 

El  número  de  multiplicaciones  y sumas  complejas  que  requiere  el 
cálculo  de  la  transformada  discreta  de  Fourier  de  la  expresión: 

F(u)  - l/N  X f<*j  cxpC-j2jBUx/N] 

i-O 


es  proporcional  iN>  que  son  necesarias  N multiplicaciones  complejas 
de  la  expresión  f(x)exp[-j2]Tux/N]  y (N-J)  sumas.  Por  otro  lado,  los  valo- 
res  de  la  exponencial  compleja  pueden  calcularse  una  sola  vez  y almace- 
narse en  una  tabla  para  posteriores  cálculos. 

Se  ha  propuesto  un  procedimiento  de  descomposición  denominado 
t musió  imada  rápida  de  fourier  (FFT),  que  reduce  el  número  de  multipli- 
caciones y sumas  a NtogjN,  lo  cual  representa  un  coste  computad  onal  mu- 
cho menor  que  el  que  requiere  la  transformación  directa*  en  especial  cuan- 
do N es  grande,  como  puede  comprobarse  en  la  siguiente  tabla 


N 

NJ 

Nlog,N 

2 

4 

2 

4 

16 

8 

S 

64 

24 

16 

256 

64 

32 

1024 

160 

64 

4096 

384 

12» 

I63S4 

896 

256 

65536 

2048 

512 

262144 

4608 
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La  mayor  aplicación  de  las  transformadas  está  en  el  diseñó  de  filtras 
para  mejora  de  imágenes  y también  en  algunas  técnicas  de  compresión  de 
imágenes.  En  el  apéndice  l se  describe  con  más  detalle  el  algoritmo  que 
permite  el  cálculo  de  la  FFT,  Asimismo,  este  algoritmo  puede  ser  utilizado 
para  evaluar  la  transformada  de  una  imagen  2-D,  efectuando,  en  primer  lu- 
gar, una  transformación  por  filas  y posteriormente  por  columnas.  En  for- 
ma de  pseudocúdigo  seria: 

ALGORITMO  FFT2D 
FOR  i:-l  TO  N-l  DO 
BEGIN 

X(i)¡-f(i,j) 

FFT(N,X) 

END; 

FORÍl-I  TON- 1 DO 
BEGIN 

FFT(N,X) 

HU  ):*X 

END. 

Cada  llamada  para  calcular  la  transformada  rápida  de  Fourier  compor- 
ta un  coste  computacional  proporcional  a Nlog,N.  Puesto  que  existen  2N 
de  tales  llamadas,  el  coste  computacional  total  de  la  FFT  bid  intensión  al  es 
proporcional  a NilogJ(NIJ. 

Este  procedimiento  supone  un  gran  avance  para  la  computación  de  la 
transformada  de  Fourier,  sobre  todo  cuando  N es  un  número  muy  grande. 
Por  ejemplu,  una  i mp  fomentación  directa  de  la  transformada  de  Fourier 
para  N - 3192,  necesita  deí  orden  de  45  minutos  en  una  máquina  como 
un  IBM  7094;  el  mismo  trabajo  realizado  con  la  implementación  de  k 
transformada  rápida  de  Fourier  emplearla  cinco  segundos  en  k misma  má- 
quina. 
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TECNICAS  DE  PROCESADO 
DE  IMAGENES 


UAS  principales  técnicas  de  procesado  de  imágenes  pueden 
dividirse  en  dos  grandes  categorías:  métodos  en  el  domi- 
nio de  la  frecuencia  y métodos  en  el  dominio  espacial.  Los 
primeros  se  basan  en  efectuar  la  transformada  de  Fourier 
de  una  imagen,  mientras  que  los  segundos  manipulan  di- 
rectamente los  pixels  de  la  imagen. 

El  fundamento  de  las  técnicas  en  el  dominio  de  la  fre- 
cuencia es  el  teorema  de  convolución.  Sea  g(x,y)  una  ima- 
~ gen  formada  por  la  convolución  de  una  imagen  f(x,y)  y un 
operador  invariante  espacialmente  Mx.y).  es  decir: 
g(*,y>  - h(x.y)  * P(x,y) 

Entonces,  por  el  teorema  de  convolución.  tenemos  la  siguiente  rela- 
ción en  el  dominio  de  la  frecuencia; 

G(u,v)  - H£u,v}  Ffu.v) 

donde  G,  H y F son  las  transformadas  de  Fourier  de  g,  h y l rcspeetiva- 
menie.  U transformada  H(ü.v)  es  a veces  llamada  función  de  transferen- 
del  proceso.  Ambas  ecuaciones  son  muy  utilizadas  en  mejora  y restau- 
ración de  imágenes,  siendo  más  impórtame  la  segunda,  ya  que  la  convo- 
I ución  discreta  es  a menudo  realizada  mis  eficientemente  en  el  dominio 
de  k frecuencia,  medíame  un  algoritmo  FFT. 

En  un  problema  típico  de  mejora,  f{x,y)  es  dada  y el  objetivo  tras 
calcular  F(u,v).  es  seleccionar  un  H(u,v)  tal  que  la  imagen  deseada,  dada 
por; 

gfx.y)  - F 1 [H(u,v)  F(utv)] 

exhiba  alguna  característica  importante  de  f{x,y).  Por  ejemplo,  tos  bordes 
en  la  imagen  f(x,y)  pueden  ser  acentuados  utilizando  una  función  H(u,v) 
que  enfatice  los  componentes  de  alta  frecuencia  de  F(u,v). 
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Las  técnicas  en  el  dominio  espacial  se  basan  en  transformaciones  de 
los  niveles  de  gris,  donde  el  tipo  de  operación  utilizada  usada  depende  dd 
criterio  empleado  para  el  realce  o mejora,  operando  siempre  sobre  el  va- 
lor de  los  pixels  de  la  imagen. 

Por  otra  parte,  podemos  hacer  otra  división  de  los  métodos  de  mejora 
de  imagen  en  el  dominio  espacial  atendiendo  al  número  de  pixels  de  la 
imagen  original,  de  los  que  depende  d valor  de  cada  pixel  de  la  imagen 
resultante: 

— Operaciones  de  punto  o pixel  individual;  En  ellas  el  valor  de  cada 
pixel  de  la  imagen  resultante  es  función  únicamente  del  valor  que  tenia 
ese  mismo  pixel  en  la  imagen  original,  mediante  una  transformación  de  al- 
gún tipo.  Dentro  de  este  grupo  se  encuentran  las  técnicas  de  modificación 
del  contraste  por  ecuatización  del  histograma  de  niveles  de  gris, 

— Operaciones  de  área;  En  ellas  se  calcula  d valor  de  cada  pixel  como 
una  función  de  los  valores  de  los  pixels  de  un  cierto  enlomo  de  dicho  pixel 
en  la  imagen  original  y no  sólo  del  valor  original  de  ese  pixel.  En  esta  ca- 
tegoría se  encuadran  las  operaciones  de  cottvolución  y algunas  transfor- 
maciones geométricas. 


MANIPULACION  DEL  CONTRASTE 

La  manipulación  del  contraste  es  una  transformación  radiometría 
puntual  (pixel  a pixel),  cuyo  objetivo  es  lograr  una  mejor  discriminación 
óc  imágenes  con  bajo  contraste.  Cada  nivel  de  gris  de  la  imagen  se  modi- 
fica mediante  la  transformación  seleccionada,  sin  considerar  los  valores 
de  los  niveles  de  gris  de  los  pixels  adyacentes.  El  análisis  de  los  histogra 
mas  de  una  imagen  constituye  en  la  mayoría  de  las  ocasiones  un  paso  pre- 
liminar para  lograr  una  eficiente  manipulación  del  contraste  de  una  ima- 
gen. El  histograma  de  una  imagen  proporciona  el  porcentaje  dd  número 
total  de  pixels  que  presentan  un  valor  determinado  de  nivel  de  gris,  El  con- 
cepto es  análogo  al  de  densidad  de  probabilidad  que  se  utiliza  en  estadís- 
tica. Las  características  generales  de  los  histugramas  de  niveles  de  gris  en 
el  caso  de  algunos  tipos  de  imágenes  aparecen  en  la  figura  1 , 


íoo  hajjt  radiancia 


A 


SI  vil  de  grlt 


Eacena  enn  illa  radiancia 


A 


Nivel  de  gri* 
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Nivel  de  gris  Nivel  de  grli 

Frjj.  /.  CanUléristica.'í  de  los  hisSo^ramdi  de  diferentes  tipos  dé  cscétitis. 


El  histograma  de  una  imagen  sólo  especifica  el  número  total  de  pixels 
de  cada  nivel  de  gris  y no  proporciona  información  acerca  de  la  distribu- 
ción espacial  de  los  niveles  de  gris.  La  figura  2 presenta  el  resultado  de  la 
aplicación  de  una  transformación  de  los  niveles  de  gris  en  el  histograma 
de  una  imagen. 


I IL|  J Resultado  dé  tu  aplicación  de  uná  transformación  fint  a!  en  el  histo^ramn. 


Cada  nivel  de  gris  de  la  imagen  de  entrada  es  proyectado  en  la  curva 
Je  transformación,  proporcionando  un  nuevo  valor  en  la  imagen  de  sali- 
da, conservándose  el  área  dd  histograma  obtenido.  Esta  simple  operación 
lineal  de  mejora  del  contraste  se  denomina  compresión  del  contraste.  En 
* p asiones,  y con  objeto  de  conseguir  un  mejor  contriste  en  la  imagen  dé 
salida,  se  sude  considerar  algún  tipo  de  saturación  en  ambos  exiremos  del 
i migo  de  salida,  como  se  muestra  en  la  figura  3, 

Si  el  histograma  es  asimétrico,  como  sucede  frecuéntenteme,  es  real- 
mente complicado  poder  controlar  mediante  una  simple  transformación 
lineal  tanto  d valor  medio  del  nivel  de  gris  de  la  imagen  de  salida  asi  como 
la  cantidad  de  saturación  en  los  extremos  de  la  transformación  efectuada 
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Fig-  t \h: jura  lineal  del  contraste  con  saturación  variable. 


Es  en  estos  casos  en  donde  podría  aplicarse  una  transformación  lineal  por 
tramos  como  Ja  indicada  en  la  figura  4 


Fifr  4.  Mejora  del  contraste  no  lineal. 


Una  operación  de  manipulación  del  contraste  más  general  lo  constitu- 
ye  Ja  denominada  ccualizacion  del  bis  (ograma.  que  se  utiliza  en  una  muy 
numerosa  cantidad  de  imágenes.  En  este  caso,  se  utiliza  como  función  de 
transformación  la  distribución  acumulativa  de  niveles  de  gris  (el  número 
de  puteis  en  el  bi  si  ograma  entre  cero  y cada  nivel  de  gris).  La  ccua- 
io.u  ión  del  hlsiograma  produce  una  reducción  det  contraste  en  áreas  muy 

J4 


clara»  o muy  oscuras,  expandiendo  los  valores  intermedios.  Un  ejemplo  de 
e<  utilización  del  hisiograma  está  representado  en  la  figura  5. 


fig.  S.  Ecuaiiiflción  del  lustograma. 


Entre  los  tipos  de  operaciones  que  pueden  incluirse  dentro  de  las  de 
manipulación  del  contraste  se  encuentra  la  denominada  como  umbral  iza- 
ciúm  El  objetivo  de  la  umb  religación  no  es  mejorar  el  contraste  de  la  ima- 
gen considerada,  sino  segmentarla  en  dos  clases  diferenciadas:  objetos  y 
futido,  mediante  la  simple  aplicación  de  un  umbral  en  el  rango  de  los  ni- 
veles de  gris.  Para  llevar  a cabo  la  operación  de  umbral ización  supónga- 
me que  los  niveles  de  gris  de  una  imagen  dada  f(x,y)  presentan  el  histo- 
grama  de  Ja  Figura  ó. 


Ftg.  6.  Vmhraliifldón  del  tas  fórrame. 

I>e  la  observación  de  la  figura  puede  verse  que  la  mayoría  de  los  pixeJs 
de  f(x,y)se  encuentran  concentrados  en  la  zona  más  oscura,  estando  el  res- 
to repartidos  más  o menos  uniformemente.  Este  histograma  es  caracterís- 
tico de  las  imágenes  que  contienen  objetos  grises  en  un  fondo  oscuro.  Para 
fia  llar  la  frontera  entre  los  objelos  y el  fondo,  dividiremos  el  hisiograma 
en  dos  zonas  separadas  por  un  umbral  T,  Se  Irata  de  seleccionar  T de  tal 
ínrroa  que  la  banda  SI  contenga  los  niveles  asociados  al  fondo,  mientras 
que  la  banda  02  contenga  los  niveles  de  los  objetos,  Asi,  un  cambio  en  el 
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nivel  de  gris  de  una  banda  a oirá  indicará  la  presencia  de  una  frontera. 
Para  conseguir  3a  detección  de  fronteras  Unto  en  sentido  horizontal  como 
vertical,  será  necesario  explorar  la  imagen  por  filas  y por  columnas.  Asi 
pues,  los  pasos  a seguir  serian  los  siguientes: 

a}  Para  cada  fila  de  Ffx.yb  crear  una  fila  de  una  imagen  intermedia 
gl(x,y),  utilizando  la  siguiente  relación  para  y»Ü,l N-L 

LE  si  los  niveles  de  gris  de  f(x,y)  y de  f(x.y-l)  glfx.y)  están  en  diferentes 
bandas 

LB  si  se  encuentran  en  la  misma  banda 

b)  Para  cada  columna  de  F(x,y),  crear  una  columna  de  una  imagen  in- 
termedia g2(x,y),  utilizando  la  siguiente  relación  para  x-0 N-J: 

LE  si  gl(xfy>  o g2(x,v)  son  iguales  a LE  g2fx,y) 

LB  si  no  son  iguales  a LE 


SUAVIZADO  DE  IMAGEN 

Las  operaciones  de  suavizado  se  usan  principalmente  para  reducir  efec- 
tos espúreos  que  pueden  presentarse  en  una  imagen  digital  como  conse- 
cuencia de  un  sistema  de  muestreo  pobre  o del  canal  de  transmisión.  Pue- 
den realizarse  tanta  en  el  dominio  espacial  como  en  el  de  Fourtcr. 

a)  En  el  dominio  espacial.  Promedio  de  un  entorno. 

Dada  una  imagen  F(x,y)  se  intenta  generar  una  imagen  suavizada  gíx,y) 
cuyos  niveles  de  gris  en  cada  pumo  fx,y)  se  obtienen  promediando  los  ni- 
veles de  gris  de  los  pixels  de  f contenidos  en  un  entorno  predefinido  de 
(Xy).  En  otras  palabras,  la  imagen  suavizada  se  obtiene  usando  la  relación: 

,m)  e S 

para  x,y  - 0.1 ,2,...rN- 1 S es  el  conjunto  de  coordenadas  de  vecinos  del  pun- 
to íx,y).  pero  sin  incluir  a éste,  y M es  el  número  total  de  estos  puntos. 

El  tamaño  del  entorno  elegido  es  arbitrario,  pero  hay  que  tener  en 
cuenta  que  el  proceso  de  suavizado  de  bordes  puede  producir  un  efecto 
de  -blumng*  o emborrona  miento  en  la  imagen,  que  no  es  deseado.  El  gra- 
do de  eriiburronamicniu  será,  además,  fuertemente  proporcional  al  radio 
di  1 emurno  usado.  Para,  un  radio  dado,  ti  efecto  de  emborronan!  iento  puc- 
di  m i i educido  usando  un  procedimiento  de  umbral ización.  Es  decir,  en 
Nqnu  de  mm  la  ecuación  anterior  para  calcular  g(x,y)  eu  todos  los  puntos, 
UMirriiHtt  dicha  expresión  si 

F(x,y)  - > T 


y en  caso  contrano  diremos  que  g(x,y)  *■  F(x,y):  donde  T es  uti  umbral  no 
negativo  predefinido.  De  esta  forma  podemos  reducir  el  nudo  de  la  ima- 
gen evitando  en  gran  parte  el  efecto  indeseado  de  ^blurring*. 

bf  En  el  dominio  de  Fouricr.  Filtros  paso  bajo, 

Los  bordes  y otras  transiciones  pronunciadas  (como  tJ  ruido)  en  los  ni- 
veles de  gris  de  una  imagen,  contribuyen  grandemente  al  contenido  de  alta 
i rv  cu  ene  iñ  de  su  transformada  de  Fourier.  Por  tanto,  el  embcuTonatniento 
pndria  ser  conseguido,  Via  dominio  de  la  frecuencia,  atenuando  un  delcr- 
n nnado  rango  de  componen! es  de  alta  frecuencia  en  la  transformada  de 
una  imagen  dada. 

Recordemos  la  ecuación: 

G(u,y)  - H(u.v)  F(u,v) 

donde  F(u,v)  es  la  transformada  de  la  imagen  que  queremos  suavizar  El 
problema  es  determinar  una  función  H(u,vj  que  nos  produzca  G(u,v>  por 
atenuación  de  las  componedles  de  alta  frecuencia  de  Ffu.v).  La  transfor- 
mada inversa  de  G(u,vj  nos  dará  posteriormente  la  imagen  suavizada  de 
sC  ada  g(x,y}.  Ya  que  las  componentes  de  alia  frecuencia  son  *el  i minadas» 
o filtradas  y la  información  en  el  rango  de  bajas  frecuencias  es  mantenida 
sin  atenuación,  este  mélodo  se  conoce  comúnmente  como  filtrado  de  pa- 
vu  bajo.  La  función  H(u,v)  se  denomina  función  de  transferencia  del  filtro 
, onsiderado.  En  todos  los  casos  que  se  considerarán  más  adelante,  las  ope- 
raciones de  filtrado  afectan  de  la  misma  manera  a la  parte  real  y a la  parte 
imaginaria  de  la  transformada  de  Fourier. 


Ejemplos  de  filtrados  paso-bajo 

Los  ejemplos  siguientes  corresponden  a filtros  paso-bajó  con  simetría 
rotacional,  es  decir,  el  filtro  completo  se  generarla  rotando  3ó0‘  la  sec- 
ción mostrada  en  cada  uno  de  Eos  ejemplos  señalados. 

a)  Rectangular,  Un  filtro  paso-bajo  cuya  función  de  transferencia  sa- 
i istaga  la  relación 

u, oí  í I si  R — RQ 
H(R)  \ o si  r > m 

Gráficamente,  en  la  figura  7 se  representa  el  perfil  de  este  tipo  de  filtro, 
que  suele  denominarse  * ideal*,  puesto  que  deja  pasar  todas  aquellas  fre- 
cuencias situadas  dentro  de  un  circulo  de  radio  menor  que  R^sin  atenuar 
las,  eliminando  lodas  aquellas  frecuencias  situadas  en  el  exterior  del 
ctrculo.  Al  valor  R^,  se  le  denomina  frecuencia  de  corte,  que  constituye  un 
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SESEir*  Ú“*  * '"hora  de  delínir  esmcteriitiCM  de  un  fll.ro,  per- 
m Hiendo  Ja  comparación  con  otros  tipos  de  filtros  ^ 

nnrn /n?oUH.d|°  ?*  k 1 *p,if“iún  * »P=> de  Blln»  produce  un  emborro- 

deZ  eb.?b¿.'ma®en-dlf  .parl'<Jí'  ,:<í,,,0  consecu«'c¡*  ¿e  I*  eliminación 

* £ Wd^™£.'0  qUe  “ '°  mlSm°'  UnB  <*«"*  d‘ 

b¡  Filtro  de  Butterworth 

La  función  de  transferencia  de  un  filtro  de  Butterworth  de  orden  n v 
de  frecuencia  de  corte  R^,  viene  dada  por  la  siguiente  expresión: 

H[R)mww^w 

Gráficamente,  en  la  fijpra  g se  presenta  el  perfil  de  este  filtro,  pudiendo 
cámo  función  de  transferencia  no  presenta  un  cambio  ian 
abrupto  como  en  el  anterior  caso  de  filtro  ideal,  siendo  en  general  mis 
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aconsejable  la  uüliiación  de  este  tipo  de  filtros  que  presentan  una  caída 
rmi*  suave,  puesto  que,  como  puede  demostrarse,  la  utilización  del  caso 
ideal  produce  falsos  contornos  en  3a  imagen  filtrada.  Seguidamente  vere- 
mos algunos  otro®  ejemplos  de  filtros  que  presentan  una  caída  suave  de  La 
irspuesla  en  frecuencias. 


ci  Filtro  trapezoidal 

tiste  tipo  de  filtro  constituye  un  compromiso  entre  el  filtro  ideal  y un 
filtro  con  caída  suave.  La  función  de  transferencia  de  este  filtro  viene  dada 
por  la  siguiente  expresión: 


si  R.  < RÜ 


H(R}« 


(RdilJ  sí  RO  S R 5 Rl 
O si  R >;  R1 


Gráficamente,  el  perfil  de  este  tipo  de  filtro  viene  representado  en  la  figu- 
ra 9, 


Fig.  9,  Filtro  paso-bajo  trapezoidal, 
df  Filtro  exponencial 

Este  tipo  de  filtro  paso-bajo,  que  también  se  utiliza  frecuentemente  en 
procesamiento  de  Imágenes,  presenta  una  frecuencia  de  curte  situada  a dis- 
tancia RO  y una  función  de  transferencia  dada  por  la  siguiente  expresión: 

H(R>-  cxp[(R/RO)n] 
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Gráficamente,  cE  perfil  de  este  lipa  de  filtros  viene  representado  en  l:i 
figura  10. 


REALCE  DE  BORDES 

Puesto  que.  como  ya  se  ha  señalado,  los  cambios  abruptos  de  los  nive- 
les  de  gris  de  una  imagen  están  asociados  con  las  componentes  de  alta  fre- 
cuencia, la  realización  de  operaciones  de  filtrado  paso-alto  producirá  una 
atenuación  de  las  componentes  de  baja  frecuencia  sin  modificar  sustan* 
l ¡símente  la  m Formación  de  frecuencias  altas  contenida  en  el  espectro  de 
Fourier.  Los  ejemplos  mis  usuales  de  filtros  paso-alio,  que  presentan  sí- 
meiría  rotacional,  son  análogos  a los  anteriormente  estudiados: 


Ejemplos  de  filtrados  paso-alto 

a)  Filtro  rectangular 


Este  tipo  de  filtros  presenta  una  función  de  transferencia  dd  tipo; 


um 


I 0 sí  R ^ R0 

I I si  R > RO 


stendo  Rij  b frecuencia  de  corte.  Gráficamente,  este  tipo  de  filtros  viene 
representado  en  la  figura  II. 


b>  Filtro  de  Buticrworth 

lut  función  de  transferencia  de  un  filtro  paso-alto  de  Butierworth  de  or- 
den n y frecuencia  de  corte  RO  viene  dada  por  la  siguiente  expresión: 


H(R)  " TT  [ kó/'R.)  Jn 
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Fi^  II.  Filtro  paso-ah®  ideal 

Gráficamente,  este  tipo  de  filtros  viene  representado  mediante  la  figu- 
ra 12. 


H|RJ 


fig.  12.  Filtro  püM^tilto  de  Bwitwu’rth. 


c)  Filtro  exponencial 

La  función  de  transferencia  de  un  filtro  paso-alto  exponencial  con  fre- 
cuencia de  corte  R0,  viene  dada  por  la  siguiente  expresión: 

H(R)“  expf-RO/R]" 

en  donde  el  parámetro  ri  controla  la  velocidad  de  subida  de  la  función  ex- 
ponencial. Una  gráfica  del  perfil  de  este  tipo  de  filtros  viene  dada  en  la  fi- 
gura 13. 
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Ptg,  13.  Filtra  pasít-aito  exponencial 
d)  Filtro  trapezoidal 

función  de  transferencia  de  un  filtro  paso-aJto  trapezoidal  viene  dada 
por  la  siguiente  expresión; 

0 si  R < R1 

H(R>  * 71CTTT  [R"R1]  *i  R1  R ^ RO 
! sí  R >-  RO 

en  donde  RO  y RJ  son  parámetros  conocidos  y se  cumple  que  R0>0R1+ 
Gráficamente,  este  tipo  de  filtros  viene  representado  en  la  figura  14. 


Hg.  / 4.  Filtra  pasQ-irf$Q  irapu  zoidat. 


32 


TRANSFORMACIONES  GEOMETRICAS 


En  ocasiones  es  necesario  efectuar  cambios  geométricos  en  una  ima- 
gen digital,  como,  por  ejemplo,  cambiar  de  sistema  de  coordenadas,  efec- 
lunr  rotaciones,  etc.  Una  veí  realizada  La  transformación  de  coordenadas 
c\  necesario  efectuar  un  proceso  de  interpolación  con  objeto  de  ajustar  la 
gen  resultante  a la  retícula  x.y.  en  donde  se  encuentra  situada  la  ima- 
gen originad-  Por  ejemplo,  la  transformación  de  coordenadas  de  un  cua- 
drilátero definido  por  cuatro  puntos  puede  efectuarse  a través  de  los  si- 
guientes polinomios; 

x-  83  + a,  *'  + / + a,*y  + a4{x')3  + a3  (y)J 

y-  bv  + bt  xl  + h2  y’  + b,x  Y + b^x'}1  + ba  (y')1 

en  donde  las  coordenadas  con  prima  son  las  de  la  imagen  de  salida  y las 
coordenadas  sin  prima  corresponden  a la  imagen  original  Si  los  coeficien- 
tes de  (xh)J  e (y1)1  son  cero,  los  ocho  restantes  coeficientes  pueden  calcu- 
larse a partir  de  los  cuatro  puntos  correspondientes  a los  vértices  del  cua- 
drilátero de  entrada.  Los  tipos  de  transformaciones  posibles  en  este  últi- 
mo caso  vierten  representadas  en  la  figura  15. 


DESPLAZAMIENTO  f~ 
* <10  + i’ 

y - bo  * f 


&ESGAMIENW 
x - x"  + a¡y' 

y m ym 


Á 


~7  PERSPECTIVA 

/ * - «?y  j 

' J 

y * y 

GIRO 

x-  *t*'  * fljy' 
y * b + ftjy' 
a,  - ttj-  Cíu  í 
ítj  - 3>t  - n 
n 


Fig.  15.  Ejemplos  de  transfonnadas  geométricas, 
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VISION  ARTIFICIAL 


A mayor  parle  dd  desarrollo  en  visión  artificial  se  ha  lie' 
vado  a cabo  en  aplicaciones  industríales  en  donde  el  en- 
lomo en  que  se  encuentra  la  máquina  es  bien,  conocido  y 
la  tarea  encomendada  al  sistema  de  visión  está  claiantcn 
te  establecida.  Un  ejemplo  típico  lo  constituye  el  sislema  di- 
visión utilizado  por  un  robot  para  manipular  objetos  en 
movimiento  situados  en  una  cinta  transportadora.  En 
otras  áreas  de  aplicación  en  donde  la  información  a ma- 
nipular está  peor  definida,  el  desarrollo  llevado  a cabo  es 
Lümparativamente  menor,  como,  por  ejemplo,  sucede  en  el  caso  del  aná- 
lisis c interpretación  de  imágenes  de  rayos  X, 

El  objetivo  de  un  sistema  de  visión  artificial  es  analizar  las  imágenes  y 
producir  descripciones  de  dichas  imágenes.  EJ  sistema  de  visión  puede  ser 
i < tusiderado  como  un  elemento  de  un  bucle  de  reali mentación  que  está  re- 
Lición ado  con  la  captación  y análisis  de  imágenes,  mientras  que  otros  ete* 
memos  deí  sistema  están  dedicados  a la  loma  e implementación  de  dichas 
decisiones  (figura  1). 


Fifi.  /.  Pane*  constituyentes  de  un  sistema  de  visión  artificial. 
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Existen  diferentes  niveles  de  descripción  posibles.  En  general,  se  con- 
siderará suficiente  aquel  nivel  de  descripción  que  permita  llevar  a cabo  la 
tarea  deseada. 

Entre  las  imágenes  más  sencillas  que  puede  tratar  un  sistema  de  visión 
se  encuentran  las  denominadas  imágenes  binarias.  Una  gran  parte  de  las 
imágenes  industriales  que  se  procesan  son  binarias,  puesto  que  en  muchas 
ocasiones  las  descripciones  que  pueden  extraerse  de  este  tipo  de  imáge- 
nes son  suficientes  para  este  tipo  de  aplicaciones. 


Una  imagen  binaria  viene  definida  por  una  función  característica 
que  es  cero  para  aquellos  valores  correspondientes  ai  fondo  de  la  imagen 
y uno  para  los  puntos  correspondientes  al  objeto.  Un  ejemplo  de  la  ima- 
gen binaria  está  representada  en  la  figura  2. 


MU) 


I 


1 


ir, 

fig-  2,  Eftmpío  de  imagen  binaria  /,  indicando  tos  centros  de  gravedad  de  cada  ab- 
je io. 


Esta  imagen  binaria  puede  proceder  de  una  operación  de  umbral  ila- 
ción como  las  anteriormente  descritas.  La  gran  ventaja  de  las  imágenes  bi- 
naria* es  que  son  fáciles  de  di  git  alizar,  almacenar  y transmitir,  puesto  que 
con  muy  pocos  bits  es  posible  representar  este  tipo  de  imágenes 

El  procesamiento  de  imágenes  binarias  requiere  tan  sólo  el  cálculo  de 
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dgunas  propiedades  geométricas  y topologías.  Veamos  seguidamente  al- 
gunas de  ellas: 

AREA:  A - 

N-1  N-l 

X X b(¡  j) 
1-0 

(4.1) 

H-f  NI 

«,-l/A  X£ib(U) 

(4.2) 

POSICION: 

y„-l/A  ££jb<ij> 

(4-3) 

Las  expresiones  de  la  coordenada  de  la  posición  corresponde  a lo  que 
se  denomina  momentos  de  primer  orden  y coincide  con  el  concepto  de 
centro  de  gravedad  de  un  objeto  de  ía  Mecánica- 

Asi  como  para  el  cálculo  de  la  posición  es  necesario  utilizar  un  mo- 
mento de  primer  orden,  para  obtener  la  orientación  es  necesario  el  cálcu- 
lo de  momentos  de  segundo  orden-  Supongamos  que  se  desea  calcular  la 
orientación  del  objeto  que  se  encuentra  situado  en  la  figura  3.  Para  calen- 


Fig.  i-  Orientación  de  un  objeto  binaría. 
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lar  d eje  según  el  cual  se  encuentra  orientada  la  figura,  hay  que  tener  en 
cuenta  que  este  eje  coincide  con  el  denominado  eje  de  mínima  inercia  que 
se  utiliza  asimismo  en  Mecánica.  Por  ello,  habré  que  minimizar  h siguiem 
te  expresión:  * 

N I NI 

E * X X r%(i  jj  (4,4) 

1-0  J-“0 


en  donde  res  b perpendicular  trazada  desde  (ij)  a b linca  anterior  (figura 
4).  Puesto  que  d eje  de  mínima  inercia  pasa  por  d centro  de  gtfvedad,  es 


F$J  Distancia  de  un  punto  US  « una  Una:  UcrA  - /«**  + 0,  ír Wa  0 la  inclina, 
ilón  respecta  ai  eje  r y 0 es  tu  perpendicular  a la  linea  desde  el  oriyen. 

conveniente  efectuar  un  cambio  del  sistema  de  coordenadas  a uno  nuevo 
que  pase  por  dicho  punto  y sea  paralelo  a!  antiguo: 
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í‘  - i-i,:  y - j-jfl 


(4.5) 


sustituyendo  en  b expresión  (4,4)  y teniendo  en  cuenta  b distancia  de  un 
punto  (L,j)  a una  recta  definida  por  ios  parámetros  & y 8.  de  b Fig,  4 se  ob- 
tiene la  siguiente  expresión: 

N-l  N-l 

E - £ £ (1‘  sen  8-  jh  eos  8 y b(Íj) 

i*0  j*0 

(4-6) 

\ desarrollando  la  expresión  anterior,  queda: 

EHa  senJ6-2b  sen(teos8+C  eosJ0 

(4.7) 

en  donde: 

i*II  b(i  j)  <>T 

1-0  r« 

(4-8) 

b - L Z b(M)  GT> 

i-tl  H9 

(4.9) 

c m Z Z b(i  j)  til1 

(4.10) 

»-0  H3 


para  hallar  La  orientación  80  dd  eje  de  mínima  inercia  habré  que  derivar 
e igualar  a cero  la  expresión  (4.7),  obteniendo  como  resultado: 

8*-  l ¡1  tang 1 (Zb/ a-c)  (411) 

En  la  Figura  5 se  presenta  un  procedimiento  en  lenguaje  USP  deno  mi 
nado  CENTRO,  que  determina  el  centro  de  gravedad  de  una  imagen  rec- 
tangular de  dimensión  N filas  x M columnas,  devolviendo  los  valores  L^,  y 
Jo  de  las  coordenadas  dd  centro  de  gravedad.  En  la  figura  6 se  presenta 
un  procedimiento  LISP  que  devuelve  una  lista  que  contiene  el  centro  de 
gravedad,  así  como  d valor  dd  ángulo  que  íorrna  el  eje  de  mínima  inercia 
con  d eje  (i). 


taenju  cehteü  t imagei 

ILET  < íH  I ARHAy-0 1 «CU5( OH  JMAtSt  > 

m iftftSftv^iMEHG  iou  i&m  Al  E 

(SUM  ai  ÍBUM-l  fKiW-J  vt 

(DO  tu  a I*  I nu 


:n  fita»* 

I»  eerlumnw** 


a*  i—  ni*c 


das  lac  fiue. 
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Fie.  5 Fu«ci<>rt  que  permite  calcular  el  área  de  une  imagen  binaria  de  dimen- 
sión jV  X W. 
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FiZ  6 Función  USP  que  permite  el  cálculo  de  la  orientación  de  una  imagen  bina 
rio  N K M 
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Otras  características  geom  Él  ricas  de  interés  pueden  ser  el  perímetro  de 
lo  figura  estudiada,  que  puede  estimarse  como  el  número  de  ceros  que  se 
encuentran  contiguos  con  unos.  En  k figura  7 se  presenta  un  procedí rruen- 
i.»  que  permite  calcular  el  perímetro  de  una  imagen  binaria,  si  bien  cabe 
« ñalar  que  el  valor  obtenido  es  una  sobreestimación  del  valor  real,  debí' 
do  □ que  el  proceso  de  cuantifkación  produce  en  la  frontera  de  los  obje- 
to* un  efecto  en  diente  de  sierra. 


{DCF1JN  PÉHtrtCTEB  Í1MASE1 

<LET  UN  AfiHftY- DIMENSIÓN  1*^  m > 

m brb ay -dimensión  iüagc  m 

<B4J*1  ÓJl 

4 00  t<  1 I < + 1 H>>  H-  I NP  SUHl 
too  tu  i í*  j un  u«  j "¡i 


in  nti»* 

|H  eelwAfi»- 
j ftctwlnj  1 «dp  ♦ 

ítln  paao 
per  nu. 

¡Un  piu  e**- 
c&lvwia . 


(COMO  4 ítwf  (-  (íwer  image  l jf 

(AfitF  niñee  i t-  j um 

4 ÜfLTO . 3uM  4*  3UH  ÍJJU 


lO-l'  t 1-0 

,'it'Uctl. 


4 CONO  UNOT  4-  tAftEF  iKflSC  I í> 

4 AflEF  I4VUBE  4-  I 1 V í >H 

4 SETO  6tm  4 + SUN  Ifl  )>»  I t 1*’“'  * 1 

FwrU^ntal. 

fíe  7 Función  USP  que  permite  calcular  el  perimetm  de  una  imagen  binaria  de 
dimensión  N X M. 


SISTEMAS  DE  PROCESAMIENTO  ^ 
DE  IMAGENES  J 


L CARACTERISTICAS  GENERALES 

OS  sistemas  de  procesamiento  de  imágenes  se  definen  como 
aquellos  capaces  de  realizar  alguna  de  las  siguientes  operar 
clon»; 

a)  Adquirir  imágenes  de  uno  o más  sensores,  represen- 
tándose éstas  en  forma  de  imagen.  La  representación  de  una 
imagen  se  efectúa  mediante  una  matriz  bi dimensional,  cons- 
tituyendo el  nivel  de  gris  de  cada  pixel  la  tercera  dimensión. 

b¡  Almacenar  imágenes  de  manera  que  puedan  ser  pro- 
cesadas, 

c)  Realizar  funciones  de  procesado  tajes  como  mejora,  restauración, 
codificación,  análisis,  etc. 

Los  componentes  necesarios  que  requieren  los  sistemas  de  procesa- 
miento para  llevar  a cabo  sus  funciones  son  los  siguientes- 

— SENSORES,  como*  por  ejemplo,  cámaras  de  TV  vid  icón,  cámaras 
de  TV  de  estado  sólido  CCD,  etc. 

_ INTERFACES  de  CAPTACION  DE  IMAGENES  (FRAME  GRAB- 

liERS), 

— MEMORIAS  DE  IMAGEN,  para  almacenamiento  y presentación, 

— PROCESADORES  DE  IMAGEN,  incluyendo  los  denominados  proce- 
sadores puntuales  y lus  procesadores  espaciales  (ariay  pruccssors), 

- MONITORES 

— MEMORIAS  DE  ALMACENAMIENTO. 

— interfaces  del  usuario  con  el  sistema. 

- SOFTWARE  BASICO,  que  permíta  al  usuario  la  realización  de  ope- 
raciones de  procesado. 


43 


Descripción  de  los  componentes  del  sistema 


a)  Sensores 

La  tecnología  estándar  para  la  imagen  electrónica  es  la  cámara  vid h- 
cón.  La  mayoría  de  3 as  cámaras  vidicón  incluyen  la  electrónica  necesaria 
para  producir  una  imagen  de  TV  convencional.  La  sensibilidad  espectral 
del  video  va  desde  350  a 700  Manómetros.  La  luz  es  enfocada  ert  un  obje 
tívo  fotosensible  que  genera  un  modelo  de  imagen  como  carga,  objetivo 
que  es  explorado  por  un  ha*  electrónico  produciendo  la  señal  de  video.  El 
ha*  electrónico  se  genera  en  un  tubo  de  vacio,  lo  cual  hace  que  este  tipo 
de  dispositivos  sean  pesados  v frágiles,  requiriendo  tensiones  elevadas 
(900  V>. 

En  i as  cámaras  que  utilizan  tecnología  de  estado  sólido  (CCD,  siglas  de 
Chargcd  Coupled  Devices),  la  luz  que  incide  sobre  cada  elemento  genera 
una  carga  proporcional  a la  intensidad  lumínica,  recogiendo  esta  carga 
mediante  condensadores,  y transferida  en  forma  de  paquetes  de  carga  a 
un  amplificador.  La  imagen  aparece  entonces  como  una  serie  de  voltajes 
a la  salida  del  amplificador.  Las  ventajas  más  importantes  de  las  cámaras 
CCD  son  su  gran  robusted  mecánica,  su  afta  sensibilidad,  difícil  saturación,, 
asi  como  un  peso  y un  tamaño  reducidos.  Sus  desventajas  más  importan- 
tes son  su  coste  elevado,  asi  como  su  baja  resolución.  En  la  anualidad  y 
debido  a sus  características,  se  tiende  a utilizar  en  entornos  industríales 
cada  vez  en  mayor  medida  las  cámaras  CCD,  puesiu  que  en  este  tipo  de  aplh 
cae  iones  la  baja  resolución  no  impide  llevar  a cabo  las  tareas  deseadas. 

b)  Interfaces  de  captación  de  imágenes  (frame  grabbersf 

En  la  mayoría  de  las  ocasiones  la  señal  de  video  procedente  de  las  cá- 
maras  de  TV  es  digitalizada  y almacenada  en  el  computador  como  una  ma- 
triz bidimensional  de  niveles  de  gris.  La  función  de  un  *frame  grabber*  es 
la  de  convertir  la  señal  de  entrada  en  un  formato  que  pueda  ser  almace- 
nado y procesado.  Este  proceso  requiere  la  digital ización  de  la  seña!  y su 
transferencia  a la  memoria  de  imagen.  Los  componentes  de  un  'ira me 
grabber*  son  los  siguientes: 

— Adquisición  de  la  señal  de  entrada. 

— Sincronización  con  el  formato  de  la  señal  de  entrada. 

” Convertidor  analógico- digital  de  alta  velocidad. 

Es  importante  tener  en  cuenta  el  formato  de  las  señales  de  entrada  a 
la  hora  de  efectuar  su  captación,  ya  que  algunas  señales,  como  las  de  TV. 
incorporan  los  sincronismos  en  la  señal  de  video,  mientras  que  en  otro 
tipo  de  señales  los  sincronismos  aparecen  exteriormentc. 
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c)  Memorias  de  imagen 

Las  imágenes  procedentes  del  sensor  se  almacenan  en  la  memoria  de 
imagen  en  posiciones  consecutivas,  representando  cada  valor  almacenado 
d nível  de  gris  de  la  imagen  considerada  Debido  al  gran  volumen  de  in 
formación  necesario  para  almacenar  imágenes,  las 

son  de  gran  capacidad  Así,  por  ejemplo,  una  imagen  de  512  x 5l2j»x*b 
x 8 bit/de  carnificación  requiere  0.25  Mbytes  de  almaccnarnien.o.  En 
cuanto  a los  tiempos  de  acceso,  es  decir,  el  tiempo  requerido  para  leer  o 
escribir  en  I*  memoria  de  imagen,  no  debe  exceder  el  requendo  para  la 
captación  de  la  señal  de  video.  Este  hecho  supone  que  en  ^nos  s sie. 
mas  de  procesado  de  imágenes  la  propia  memoria  de  refresco  se  utilice 
como  memoria  de  imagen.  . 

La  memoria  de  refresco  constituye  la  parte  más 
de  visualización,  ya  que  su  diseño  es  critico.  En  los  v.suaiizadores  remen 
el  hecho  de  que  la  memoria  de  refresco  se  utilice  como  memoria  de 
imagen,  permite  realizar  funciones  tales  como  ampliación, 

,o  corte  etc.  Estas  operaciones  no  destruyen  los  dalos  contenidos  en 
memoria  En  la  mayoría  de  los  sistemas,  la  memona  de  refresco  se  em 
en  forma  de  canales.  Un  canal  típico  contiene  una  una- 
“slÍ*  512  x 8 bits,  lo  cual  equivale  a refrescar  una 
croma  en  pantalla.  Si  se  utilizan  tres  canales  en  paralelo,  podrían  repre 
sentarse  imágenes  en  color.  Los  sistemas  actuales  proporcionan  entre  1 y 
*,4  canales,  existiendo  asimismo  sistemas  que  proporcionan  imágenes  de 
1024  x 1024  pixels  y ló  bits  de  precisión. 

El  computador  central  accede  a la  memoria  de  refresco  mediante  al- 
guno de  los  siguientes  métodos:  l)  a través  de  una  mterface  que 
acceder  al  usuario  de  forma  aleatoria  por  medio  de  una  transferencia  D . 
^ considerando  una  subregión  rectangular  del  canal 
de  direccionamiento  del  computador,  por  lo  que  se  accede  a la  m.sma 
como  si  se  tratase  de  la  memoria  principal  del  computador. 

Desde  mediados  de  los  setenta,  las  memorias  de  refresco  se  han  cons- 
truido utilizando  memorias  RAM  Las  sucesivas  mejoras  en  la  tecnólogja 
RAM  ha  permitido  incrementar  la  capacidad  de  los  chips,  por  ejemplo, 
Kbíts  16  Klbiís,  64  Kbils,  etc.  Sin  embargo,  existe  una  relación  inversa  en- 
tre el  tiempo  de  acceso  a este  tipo  de  memorias  y el  consumo  de  potencia, 
porcentaje  de  errores  y coste.  Es  por  ello  que  las  memorias  de  refresco  se 
han  construido  a partir  de  dispositivos  relativamente  lentos.  Un  valor  típi- 
co del  tiempo  de  ciclo  de  este  tipo  de  memorias  es  el  de  400  nseg.,  que  es 
mucho  más  lento  que  el  requerido  para  el  refresco  de  los  visualizadores. 
Para  obtener  las  velocidades  de  refresco  necesarias,  las  memorias  de  los 
sistemas  utilizan  multipl exación  de  pixels,  típicamente  con  un  factor  de  s 
en  los  visualizadores  que  utilizan  un  formato  de  512  * 512.  En  general,  la 
mayoría  de  los  sistemas  multiplexan  cada  bit  de  intensidad  en  un  chtp  di- 
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ferenie;  por  ello,  en  un  diseño  típico  se  accede  a 64  chips  de  RAM  durante  ] 
un  simple  ciclo  de  memoria. 

I-a  mayoría  de  las  I imitaciones  de  tas  memorias  de  refresco  proceden  ] 
de  Ea  utilización  ¡Je  la  multiplicación.  Esie  hecho  supone  que  algunos  vi-  1 
suallzadorcs  no  dispongan  de  todos  los  valores  enteros  de  ampliación  i 
tsuom),  siendo  los  valores  2:1,  4: 1 y 8:1  mucho  más  fáciles  de  implemcm  I 
tar,  La  multiplexación  también  limita  el  acceso  por  columnas,  puesto  que  ] 
los  ocho  pixels  muliiplcxados  son  direccíonados  por  filas  La  consecüen*  ] 
cia  de  dio  es  que  el  acceso  por  columnas  será  ocho  veces  más  lento  que 
por  tilas. 

d)  Monitores  de  preterí  rae  ion  de  ánágenex 

El  monitor  de  presentación  de  imágenes  constituye  un  monitor  indus-  | 
trial  de  alta  calidad.  Generalmente,  se  encuentra  calibrado  para  presentar  ] 
una  trama  cuadrada  de  datos  con  una  relación  de  aspecto  1:1  Los  tama-  I 
ños  de  presentación  estándar  son  512  x 512  o 1024  x lü24  puteis,  sí  hien  I 
también  suelen  encontrarse  otros  tamaños  tales  como  640  k 480,  I 
1 024  x 1 280,  La  señal  de  TV  comercial  utiliza  velocidades  de  30  (o  25)  ira-  I 
mas/seg..  con  una  relación  de  aspecto  2:1  y entrelazamiento.  Sin  embar-  ] 
go,  se  vienen  utilizando  recientemente  visual izadorcs  con  un  refresco  de  ] 
60  (o  50)  tramas/scg,  y sin  entrelazamiento.  Este  último  tipo  de  visualiza-  1 
dores  produce  una  presentación  carente  de  parpadeos,  lo  cual  minimiza  { 
la  fatiga  del  observador. 

Por  otro  lado,  es  bastante  frecuente  la  utilización  por  parte  de  algunos  ] 

la bríc antes  de  t erm i nologia  a mbigua h fu ndame nta! mente  en  e I caso  del  tér- 
mino resolución.  La  única  medida  de  resolución  aceptable  consiste  en  co- 
locar un  test  de  prueba  en  el  monitor  y medir  ópticamente  la  resolución 
con  un  fotómetro. 

Las  principales  restricciones  de  los  visualizado  res  se  deben  a la  diver- 
gencia dd  haz  de  electrones  y a la  dispersión  de  los  fotones.  En  la  figura 
I se  presenta  la  arquitectura  típica  de  un  visual ízador  de  un  sistema  de  pro- 
cesamiento de  imágenes- 

c)  Procesadores  de  ¡mugen 

Existen  un  gran  número  de  arquitecturas  para  procesamiento  de  imá-  j 
genes,  incluyendo  las  arquitecturas  especializadas  En  general,  cabe  seña- 
lar que  tos  procesadores  serie  convencionales  tipo  vori  Ncuman  no  son  Jos 
más  adecuados  para  d procesado  de  imágenes,  sino  más  bien  arquitectu- 
ras de  tipo  paralelo. 

Las  principales  arquitecturas  para  procesamiento  de  imágenes  son  las 
siguientes: 

— Procesadores  serie  convencionales. 

- Procesadores  con  paralelismo  lineal  (pipdine). 
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— Procesadores  paralelo  (array  proccssois). 

— Procesadores  de  propósito  especial 

— Procesadores  VLSI. 

Como  anteriormente  se  indicó,  los  procesadores  serie  no  son  los  más 
adecuados  para  el  procesado  de  imágenes.  Sin  embargo,  presentan  venta- 
jas. sustanciales  con  el  rcsló  de  procesadores,  como  son  su  bajo  cosle  y su 
tlexibilidad.- 


= Métodos  de  iluminación 

Un  aspecto  importante  a tener  en  cuenta  a ia  hora  de  efectuar  lá  cap- 
tación de  una  imagen  pur  parte  de  un  sistema  de  visión  Lo  constituye  el  mé- 
todo empleado  para  la  captación  de  dicha  imagen.  Cuando  lo  que  se  desea 
es  captar  una  escena  industrial  es  necesario,  en  primer  lugar,  decidir  el  nu- 
mero y l.i  configuración  de  tas  cámaras  a emplear,  puesto  que  si  se  desea 
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obtener  información  tridimensional,  será  necesario  el  concurso  de  varias 
cántaras.  Otro  de  los  aspectos  a tener  en  cuenta  es  el  método  de  ilumina- 
ción empleado,  puesto  que  una  iluminación  inadecuada  puede  producir 
sombran  brillos  u otros  efectos  indeseados.  Los  métodos  más  usuales  de 
iluminación  de  escenas  son  los  siguientes: 

a)  iluminación  difusa.  Es  un  tipo  de  iluminación  indirecta  que  condu- 
ce a una  buena  eliminación  de  sombras,  permitiendo  una  modclización  tri- 
dimensional del  objeto  considerado  ut Mirando  una  sola  cámara  (figura  2). 


Fig.  2.  E¡empía  de  disposición  de  una  iluminación  difusa. 


b)  j Iluminación  posterior , Constituye  el  sistema  de  iluminación  más  uti- 
lizado, Permite  el  reconocimiento  de  piezas  pequeñas  y ligeras,  por  lo  que 
es  un  método  muy  indicado  para  aplicaciones  de  control  de  calidad  y cla- 
sificación de  piezas  (figura  3), 


Fin  i. 


Disposición  de  iluminación  posterior. 


4 H 


c )  Iluminación  modulada.  Este  tipo  de  iluminación  permite  la  extrac- 
ción de  información  de  curvaturas,  asi  como  de  su  uniformidad  superfi- 
cial por  medio  de  la  deformación  que  experimentan  unas  lincas  proyecta- 
dos sobre  el  objeto,  mediante  la  utilización  de  una  máscara  conocida  (figu- 
ra 4). 


Fig  4,  Disposición  de  iluminación  modulada. 


d)  Iluminación  direccional.  Este  tipo  de  iluminación  se  utiliza  funda- 
menta I mente  para  la  observación  de  las  rugosidades  superficiales  de  un  ob- 
jetó, asi  como  para  la  observación  de  pequeñas  áreas  del  objeto  conside- 
rado* mediante  la  utilización  de  un  haz  coherente  (láser)  o no  {figura  S). 


Fig.  5 . Disposición  de  iluminación  direCcional. 
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CARACTERISTICAS  HARDWARE  Y SOFTWARE 
DE  ESTOS  SISTEMAS 


En  los  últimos  diez  años,  los  sistemas  para  procesamiento  de  imágenes 
han  experimentado  un  considerable  desarrollo,  de  tal  manera,  que  upen* 
dones  de  procesado  que  en  tos  primeros  sistemas  suponían  un  elevado 
tiempo  de  cálculo,  en  la  actualidad  se  realizan  prácticamente  en  tiempo 
real.  A partir  de  los  años  setenta,  estos  sistemas  se  han  venido  utilizando 
fundamentalmente  en  investigación-  La  razón  fundamental  del  inierés  por 
el  diseño  de  este  upo  de  sistemas  radica  en  el  hecho  de  que  la  lee  no  logia 
serie  de  los  computadores  existentes  no  es  la  más  adecuada  para  el  proce- 
samiento de  este  tipo  de  información. 

Los  primeros  intentos  de  construcción  de  este  tipo  de  sistemas  se  efec- 
tuaron en  la  década  de  los  cincuenta.  Estos  sistemas  se  caracterizaban  por 
su  gran  tamaño,  memoria  limitada  y mala  calidad  de  imagen,  debida  ai  rui' 
do.  El  procesamiento  digital  de  imágenes  comenzó  realmente  a fines  de 
los  sesenta,  con  el  desarrollo  de  computadores  y m i nicom pujadores  más 
baratos  que  los  anteriores.  Un  ejemplo  del  software  de  procesamiento  de 
imágenes  desarrollado  en  esa  época  es  el  VICAR,  construido  por  el  Jet  Pro- 
pulsión Laboratorv,  Este  tipo  de  sistemas  estaban  basados  en  un  ordena- 
dor central  y las  tareas  se  atendían  mediante  una  cola  de  trabajo  y eran  de 
un  coste  excesivamente  caro  como  para  poder  ser  comercializados.  Algo 
más  tarde,  el  Laboratorio  de  Procesamiento  de  Imágenes,  de  la  Universi- 
dad del  Sur  de  California,  efectuaba  los  primeros  procesamientos  locales 
con  ayuda  de  un  microcomputador. 

La  arquitectura  propuesta  a principios  de  los  setenta,  basada  en  un  mi- 
nicomputador junto  a una  serie  de  periféricos,  se  ha  mantenido  hasta  la  ac- 
tualidad- Los  periféricos  utilizados  permiten  la  adquisición  y d igUal  ización 
de  las  imágenes,  su  presentación  y registro.  Por  otro  lado,  y con  objeto  de 
aumentar  la  velocidad  de  procesamiento,  se  les  suele  conectar  aun  com- 
putador central  o incluso  añadirles  un  procesador  de  array. 

Los  primeros  visual  izado  res  comerciales  fueron  el  Acrojet  (posterior- 
mente denominado  Cornial)  y permitían  la  realización  de  manipulaciones 
radiomét ricas,  superposición  de  gráficos  y acceso  directo  a los  pixels  de 
b imagen  mediante  cursor.  Estos  visual  izado  res  estaban  limitados  a la  pre- 
sentación de  las  imágenes  y a manipulaciones  globales  de  la  intensidad  o 
del  colar, 

En  la  mitad  de  los  setenta,  se  desarrollaron  los  pn meros  sistemas  de 
procesamiento  realmente  interactivos  con  la  incorporación  de  software  de 
tiempo  real.  Asi,  por  ejemplo,  cabe  citar  el  sistema  VICAR,  que  posterior- 
mente  dio  lugar  al  Mini-VICAR,  el  sistema  GE  100,  de  Genera!  Electric 
pum  procesamiento  interactivo  de  imágenes  LANDSAT,  o el  sistema  FE- 
i OS,  de  Eleciromagnetic  Systems  Laboratory.  La  evolución  de  las  nuevas 


no 


i reno  logias  de  computadores  fue  acompañada  de  una  progresiva  evolu- 
ción de  los  sistemas  de  presentación.  Asi,  3a  incorporación  de  memorias 
de  acceso  aleatorio  RAM  como  memorias  de  refresco  permitió  la  realiza- 
ción de  operaciones  más  complejas,  como,  por  ejemplo,  la  suma  y resta 
efectuadas  a velocidad  de  video;  el  cálculo  de  hisi agramas  por  hardware, 
particiones  y ampliaciones  de  la  imagen,  etc. 

A mediados  de  Eos  setenta.  Intermitió (tid  Imaging  Systems  introdujo  el 
denominado  Modelo  70,  que  incorporaba  un  bucle  de  realimcniáción  a ve- 
locidad de  video,  que  con  el  paso  de  los  años  se  convirtieron  en  una  ca- 
racterística común  a lodos  los  equipos,  incorporando  unidades  animé  ti- 
fo-lógicas  para  aumentar  la  velocidad  (como,  por  ejemplo,  los  sistemas 
Gemid  ÜeAnza,  Vicom,  Ikonas,  etc).  El  procesamiento  en  forma  de  bucle 
de  rcaiimentación  refleja  una  de  las  características  básicas  de  los  algorit- 
mos de  procesamiento  de  imágenes:  h realización  de  operaciones  simples 
i un  imágenes  puede  ser  aplicada  i lenitiva  mente  a cada  uno  de  los  pixels. 

Los  vi suali ¿adores  desarrollados  a fines  de  los  setenta  y principios  de 
los  ochenta  incorporaban  las  siguientes  características:  I > operaciones  op- 
cionales para  la  realización  de  operaciones  específicas,  tales  como  opera- 
nones  de  convolución,  y 2)  arquitecturas  paralelas  y bucles  de  realimeti- 
t ación  para  la  realización  de  algoritmos  iterativos. 


Sistemas  comerciales  de  procesamiento 
digital  de  imágenes 

Los  sistemas  comerciales  de  procesamiento  digital  de  imágenes  pue- 
den ser  de  tres  tipos: 

* Sütenws  de  propósito  genera!.  Están  constituidos  por  un  minicompu- 
tndor,  junto  con  los  periféricos  necesarios  para  efectuar  las  operaciones 
de  procesamiento.  En  otras  ocasiones  se  comercializa  sólo  el  procesador 
de  imágenes  para  ser  conectado  a algún  computador  central-  Ejemplos  de 
este  tipo  de  sistemas  lo  constituyen  d sistema  VICOM  y el  GQULD  DcAnza. 
En  la  figura  ó aparece  un  esquema  de  la  arquitectura  del  sistema  de  pro- 
cesamiento de  imágenes  VICOM,  Como  caso  particular  de  este  tipo  de  sisr 
lentas  se  encuentran  los  denominados  de  propósito  específico.  Constitu- 
yen sistemas  cerrados  utilizados  generalmente  para  aplicaciones  concre- 
tas, Así,  por  ejemplo,  cabe  citar  el  sistema  InterScan,  de  Control  Automa- 
ción. En  Ja  figura  7 aparece  un  cuadro  comparalivo  de  algunos  de  los  sis- 
temas comerciales  existentes,  tanto  de  propósito  general  como  de  propó- 
sito específico. 

* Esiacionex  de  ingeniería  ( warkstútions j,  Eíte  tipo  de  sistemas  propor- 
cionan una  potencia  local  de  cálculo  elevada,  de  tal  manera  que  et  usua- 
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rio  puede  llevar  a cabo  tareas  complejas  a un  coste  moderado,  sin  preci- 
mi  el  concurso  de  un  computador  central,  Sin  embargo,  las  estaciones  de- 
dicadas al  procesamiento  de  imágenes  no  suelen  aparecer  autónomas 
í'4áfíd~akmt),  sino  conectadas  por  medio  de  una  red  local  como  la  Ether- 
net. De  este  modo,  los  usuarios  pueden  compartir  recursos  comunes  a to- 
dos los  usuarios  de  la  red,  como  son  los  dispositivos  de  almacenamiento 
dr  alta  capacidad  (filé  serven},  asi  como  dispositivos  de  salida  comunes. 
Ejemplos  de  este  tipo  de  estaciones  lo  constituyen  la  estación  SUN  y APO- 
LLO. 


* Computadores  personales.  Los  sistemas  de  procesamiento  de  imáge- 
nes basados  en  computadores  personales  son  adecuados  en  aquellas  apli- 
caciones en  donde  el  bajo  coste  de  la  instalación  sea  un  punto  fundamen- 
tal Presentan  la  gran  ventaja  de  disponer  de  un  número  muy  amplio  de 
[tu jetas  de  adquisición  de  imágenes  y de  software  asociado  Como  ejem- 
plos de  este  tipo  de  tarjetas  cabe  citar  las  de  DATA  TRANSLATlON  y las 


de  PC  VISION. 


Características  de  los  algoritmos  para  procesamiento  de 
imágenes 

Los  algoritmos  para  procesamiento  de  imágenes  deben  presentar  una 
característica  común  a todos  ellos,  como  es  su  naturaleza  no  excesivamen- 
te compleja,  puesto  que  deben  ser  aplicados  a millones  de  datos.  Como  an- 
teriormente se  indicó,  los  algoritmos  para  procesamiento  de  imágenes  pue- 
den ser  divididos  cu  dos  grandes  grupos  (figura  8): 

* Algoritmos  de  procesamiento  puntual  En  este  tipo  de  algoritmos,  las 
operaciones  se  llevan  a cabo  sobre  cada  pixcl  individualizado,  sin  consi- 
tk-rar  el  valor  de  los  pixets  próximos.  Como  ejemplos  cabe  citar  las  ope- 
raciones aritméticas  y booleanas  entre  imágenes,  las  transformaciones  ra- 
il iomét  ricas  y las  operaciones  estadísticas.  Las  iransformac  iones  radiomé- 
I ricas  modifican  tos  niveles  de  gris  de  una  imagen  por  medio  de  la  aplica- 
non  de  funciones  lineales  y no  lineales  (incluyendo  fundones  aleatorias), 
dentro  del  rango  de  precisión  del  hardware  disponible. 

Constituyen  la  base  de  una  gran  parte  de  las  operaciones  de  procesa- 
miento, incluyendo  la  modificación  del  contraste.  Las  operaciones  radio- 
métricas  se  llevan  & cabo  utilizando  las  denominadas  tablas  de  consulta 
(look-up  rabies:  LUT),  que  se  encuentran  simadas  entre  la  memoria  de  re- 
úesco  y el  monitor.  Las  tablas  de  consulta  se  encuentran  implemeniadas 
en  memoria  RAM  de  acceso  rápido  (menor  de  100  nseg.),  empleándose  la 
intensidad  de  ios  pixels  de  la  memoria  de  refresco  como  una  dirección  en 
la  tabla,  proporcionando  un  valor  de  salida  que  se  presenta  en  el  monitor. 
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fo)  Transformación  puntual 


Este  proceso  se  realiza  a frecuencia  de  video,  permitiendo  que  los  valores 
almacenados  en  la  memoria  de  refresco  puedan  ser  transformados  en  otros 
valores  y representados  en  la  pantalla  de  acuerdo  con  la  función  que  se  en- 
cuentre cardada  en  la  tabla,  AslT  por  ejemplo,  la  escala  de  gris  puede  ser 
invenida,  o efectuar  una  transformación  radiométrica  logarítmica  con  ob- 
jeto de  aumentar  d contraste  de  las  zonas. 

* Algoritmos  de  procesamiento  espacial  En  este  tipo  de  algoritmos  la 
imagen  resultante  de  las  operaciones  llevadas  a cabo  es  función  tanto  de 
las  Intensidades  de  Jos  piséis  de  entrada  como  de  la  de  sus  vecinos. 

Los  algoritmos  de  procesamiento  puntual  son  de  fácil  implementación 
en  los  computadores  convencí  un  a Jes,  mediante  la  utilización  de  operacio- 
nes de  entrada/ sal  ida,  Sin  embargo,  los  algoritmos  de  procesamiento  es- 
pacial requieren  un  mayor  grado  de  complejidad  que  los  anteriormente  se- 
ñalados. Como  ejemplos  de  operaciones  de  procesamiento  espacial  se  en- 
cuentran las  operaciones  de  convofución,  en  donde  el  resultado  se  obtie 
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Proyección  de  una  ventano  sobre  una  imagen  que  está  siendo  procesada. 


Pínol  1 ptul  díntru  de  OH*  linea  Línta  M lie*» 

{b}  Cálculo  de  ventanas  consecutivas. 

Fig.  9.  Ejemplos  de  operaciones  de  canvotución. 

ne  multiplicando  cada  pixel  de  la  imagen  original  por  un  núcleo  o venta- 
na de  N x H pixels  (figura  9).  En  la  practica,  se  utilizan  núcleos  de  3x3 
o de  5 x 5 en  las  convol liciones  en  el  dominio  espacial,  con  objeto  de  que 
et  tiempo  necesario  para  llevar  a cabo  esie  tipo  de  operaciones  no  sea  ex- 
cesivamente grande. 

Por  otro  lado,  y por  razones  de  coste,  las  operaciones  de  procesamien- 
to espacial  se  realizan  de  forma  iterativa,  procesando  un  punto  del  núcleo 
ti  ventana  por  trama,  con  lo  que  se  consigue  un  tiempo  de  ejecución  próxi- 
mo al  tiempo  real. 
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Un  aspecto  importóle  □ tener  en  cuenta  a h hora  de  evaluar  los  al* 
gori  Irnos  para  procesamiento  de  imágenes  es  el  de  su  tiempo  de  ejecución 
Podemos  considerar  tres  tipos  de  intervalos  témpora  Jes: 

a)  Tiempo  real.  Se  considera  que  un  algoritmo  para  procesamiento  de 
imágenes  se  ejecuta  en  tiempo  real  si  dicho  tiempo  es  comparable  al  ne- 
cesario para  efectuar  el  refresco  de  i as  memorias  de  imagen,  de  tal  mane- 

t^ue1e3  ^ario  apenas  aprecie  que  ha  transcurrido  un  tiempo  desde  que 
lanzó  Ea  orden  de  ejecución  y el  mu  mentó  en  que  se  presenta  el  resultado. 

b)  Próximo  al  tiempo  real.  Se  considera  que  un  algoritmo  para  proce- 
samiento de  imágenes  opera  en  un  tiempo  próximo  al  tiempo  real  si  el 
tiempo  de  ejecución  es  del  orden  de  3 segundos,  Como  ejemplo  de  este 
tipo  de  algoritmos  cabe  señalar  las  eonvoJuc iones  espaciales. 

c)  No  tiempo  real  Se  considera  que  un  algoritmo  para  procesamiento 
e imágenes  opera  en  un  tiempo  muy  superior  al  tiempo  real  si  el  tiempo 

de  ejecución  está  por  encima  de  los  3 segundos  anteriormente  señalados. 
Esto  tipo  de  algoritmos  que  suponen  un  tiempo  elevado  de  ejecución  se 
suelen  introducir  en  una  cola  de  procesos. 


SISTEMAS  DE  PROCESAMIENTO  DE  IMAGENES 
BASADOS  EN  UN  COMPUTADOR  PERSONAL 

Durante  los  últimos  cinco  afios,  los  sistemas  de  procesamiento  de  imáge- 
nes  basados  en  ordenadores  personales  han  comenzado  a jugar  un  papel  im- 
portante en  algunas  aneas  de  aplicación,  habiendo  experimentado  durante 
los  uíiimos  dos  artos  un  rápido  crecimiento  la  tecnología  de  este  tipo  de  sis- 
temas. El  progresivo  descenso  de  los  costes  tanto  del  hardware  como  del  soft- 
ware convierten  al  ordenador  personal  en  una  herramienta  muy  adecuada 
cuando  el  coste  reducido  del  sistema  es  un  requerimiento  básico, 

ua  configuración  típica  de  un  sistema  de  procesamiento  de  imágenes 
basado  en  un  ordenador  pegona!  IBM  PC-AT  con  512  Kbytes  de  memoria 
principal.  20  Mbytes  de  memoria  en  disco,  requeriría  el  siguiente  hardwa- 
re adicional: 

— Tarjeta  de  adquisición  de  imágenes. 

— Cámara  de  video  y unidad  de  visual  i ración, 

**  Un  ratón  conectado  vía  serie. 

Opcional  mente  puede  añadirse  el  siguiente  hardware; 
genes  COpn>CC5adoni  conectada  a la  tarjeta  de  adquisición  de  imá- 

Tarjeta  adoptadora  aun  bus  externo,  muhtbús,  ATr  etc. 

— Tarjeta  de  conexión  a una  red  local:  Ethernet,  PC-NET,  etc. 

— Tarjeta  de  interface  a un  disco  óptico,  etc. 
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Fi«  W. 


Estructura  de  ua  sisu™  tj*  procesamiento  de  imágenes  basado  en  un 
deitudur  personal 


5V 


Yisufthia 


de  wi“ TZ  P1™1***  «lectura  de  un  sistema  de  tratamiento 
de  imágenes  basado  en  un  ordenador  peraonal  En  cuanio  al  software  ne 
STriíft  . señalar  que  ti  mayor  inconveniente  que  presenta  este  tipo 
de^onfipyrsciones  es  b limitada  capacidad  de  direecionamiento  de  me- 
mona  de  Ja  versión  de  [ sistema  operativo  MS-DOS  3.2  Este  problema  pue- 

ve™oní™dMSUM“  OPera!ÍV°  XEN‘X  ° “ fu'1™ 

Una  de  las  mayores  ventajas  de  la  zación  de  este  tipo  de  s.slemas 

reside  en  d hecho  de  h gran  cantidad  de  hardware  y de  software  dispon!- 
Oles  que  tac  i litan  la  implcmentación  de  un  sistema  de  tratamiento  de  imá- 
genes para  un  coste  reducido. 
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APENDICE 


ALGORITMO  DE  LA  TRANSFORMADA  RAPIDA 
DE  FOURIER 

Consideremos  en  primer  lugar  la  expresión  de  la  TDF: 

ni 

F(u)-  I/N  Z f(x)  exp[-j2ftux/IM]  - 1/N  Z fí*>  (A.1) 

siendo  - exPí-j2n/M]  ,N-2nyN-2M.  Seguidamente,  introducire- 
mos un  nuevo  indice  K en  el  sumatorio,  relacionado  con  x medíanle  las 
siguientes  relacionen 

x - 2K  para  x par 
x - 2K  + l para  x impar 

con  lo  que  La  expresión  (A  l)  puede  escribirse  como: 

F(u)  - 1/2M  £ fitk>  W*  - I/2H/M  £ f(2K)  W,“®  + 

l/M  £ f(2K  + I)  w£*''l  - l/2ftyw  + F1W1(U)  W,;i  (A.2) 

En  definitiva,  se  ha  reducido  d cálculo  de  una  transformada  de  N pun- 
tos a dos  transid  miadas  de  (N/2)  puntos.  Asi,  por  ejemplo,  en  el  caso  de 
una  función  con  8 punios  lf<0),  ft  I ),  ....  f(7)f,  habría  que  considerar  por  un 
lado  Los  muestras  pares  lf(0),  f(2).  F{4),  f(6)S  y por  otro  las  impares  lf(  1).  í(3)> 
ffS),  f{7)|,  Cada  una  de  las  transformadas  de  4 punios  se  calcula  como  dos 
transformadas  de  2 puntos,  Asi.  para  calcular  la  FFT  de  11(0),  f(2),  f{4),  Í(fc)í 
hav  que  dividir  este  conjunta  en  dos:  la  parte  par  íf(Q),  f(4)l  y la  parte  im- 
par lf(2),  f(6)|.  Análogamente  se  haría  en  d casa  del  segundo  conjunto  de 
las  cuatro  muestras  impares:  |f(l),  f(3)h  f(5),  f(7)l,  que  se  dividiria  en  tf(l), 
f(5)t  y [f(3)h  í(7)h  Asi,  pues,  el  algoritmo  operarla  según  el  esquema  de  la 
figura  I „ dé  arriba  a abajo: 

Seguidamente  se  presenta  una  subruiina  FORTRAN,  propuesta  por 
Cooley  y Tuckey  en  19P9  que  permile  calcular  la  FFT  mediante  el  método 
descrita  anteriormente.  Los  parámetros  de  la  subrutina  son  dos.  F es  un 
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TRANFORMADAS  DE 
ocim  pontos 


F'e  ‘ SrpF? \CiÚ"  * ,0S  da">‘  de  la  lur“:Um  IfWli  #‘> f<7)\  pon  ti  cálculo  de 


¡L™ía  °rr"aíliJ'’c  d'*“  cakular  y LN  - n,  siendo  N - 2n  ¡a  d¡- 

“"“'‘‘"do.  Con.0  resultado  de  ¡os  cálculos  esta  sub- 
oue  F ha  nen  F e Vj0r  d'  'ransformada  de  Fourier.  Obsérvese 
™™  i ■ ? ^'"sionado  ™ 'a  subrutina  como  uti  army  de  números 

ontplejos,  luego  si  partimos  de  un  army  f correspondiente  a una  función 
' CCr0anICSde  “amara  esta subrutino. 

COM  PLEX  F(  1 024),  UF  WrTFCM  PLX 

PI-3.I4Í593 

N-2«LN 

MV2-N/3 

NMUN'l 

J-l 

DO  3 l-lPNM¡ 

[F(tCEJ)  GO  TO  I 
T*F(J) 

F(J)-FC1) 

ra-r 

1 K-MV2 

l IF(K.GE.J)  GO  TO  3 
I-J-K 


K-K/2 
GOTO  2 
Í J-J+K 

DO  5 L-I.LN 
LE“2"*L 
LEI -LE/2 
ü-íl. 0*0,0) 

W- CMPLX  ( COS(  P l/LE  I),-S1N(PI/LEI)) 
DOS  J-l.LEI 
DO  4 I-J.N.LE 
IP-HLEI 

T-F(tP>n) 

F(1P)-F(I)+T 

4 F(])-F(l)+T 

5 U-li*W 
DO  6 I=1FN 

6 F(  1 ) » F(  I )/FLO  AI  ( N) 

RETURN 

EMD 
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PROGRAMAS  DE  EXPERIMENTACION 

El  objetivo  de  los  programas  que  se  presentan  en  este  apéndice  consis- 
te en  proporcionar  una  herramienta  sencilla,  flexible  y Fácil  de  utilizar  por 
aquellos  usuarios  de  ordenadores  personales.  Los  programas  se  han  escri- 
to en  Fortran  77.  requiriendo  para  su  utilización  una  impresora  de  carac- 
teres con  objeto  de  poder  imprimir  los  resultados. 

Las  imágenes  utilizadas  tienen  un  tamaño  de  64  x 64  ptxek,  almace- 
nando cada  pixel  como  un  valor  complejo  en  punto  flotante.  Este  tamaño 
de  imágenes  es  muy  apropiado  para  su  presentación  tanto  en  un  monitor 
alfanu mérito  como  para  su  impresión  en  una  impresora  de  líneas,  propor- 
cionando suficiente  resolución  para  el  fin  didáctico  a que  se  dirige.  Este 
formato  de  datos  complejos  se  ha  efectuado  con  objeto  de  Facilitar  el  aná- 
lisis de  Fourier.  La  precisión  de  punto  flotante  asegura  la  no  introducción 
de  errores  de  redondeo.  Por  otro  lado,  el  cálculo  de  la  transformada  de 
Fourier  se  ha  efectuado  de  tal  manera  que  la  frecuencia  cero  corresponde 
al  centro  de  la  imagen  presentada.  La  presentación  en.  la  pantalla  se  ha  he- 
cho utilizando  una  sola  línea  de  caracteres  para  producir  ocho  niveles  de 
seudogris 

En  el  programa  se  han  definido  tres  buffets  de  imágenes  de  64  x 64 
pixels,  lo  cual  permite  la  realización  de  operaciones  de  suma,  escalado, 
transformación,  etc.  La  utilización  de  este  programa  de  entrenamiento  es 
sencilla,  puesto  que  se  basa  en  la  elección  entre  una  lista  de  comandos  pre- 
definidos acompañados  de  uno  o dos  argumentos  correspondientes  a los 
*buffcrs>  de  imagen.  El  resto  de  los  parámetros  se  introducen  en  respues- 
ta a las  peticiones  efectuadas  por  el  programa,  asignándose  valores  por  de- 
fecto en  el  caso  de  no  especificar  dichos  parámetros. 

Asimismo,  el  programa  permite  el  diagnóstico  de  errores  de  sintaxis  en 
los  comandos,  devolviendo  control  al  usuario  para  su  nueva  introducción. 
Los  errores  producidos  como  consecuencia  de  la  introducción  de  un  pa- 
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ramciro  erróneo  implican  la  necesidad  de  introducir  de  nuevo  el  coman- 
do, En  la  figura  I se  presenta  una  imagen  de  barras  y en  la  2 el  módulo  de 
su  transformada  de  Founer.  En  la  3 se  presenta  una  imagen  constituida 
por  un  cuadrado  de  niveles  de  gris  y en  la  4 el  módulo  de  su  iransformada 
de  Founer, 


PflOOftjin  AUTOMATA 

c 

e E*t*  pr(b?r**«  #ft  *1  qu«  p0r  WfflI> 

C de  AutümM*  L*i  op*r*eÍDnc*  fluí  „ 

C p^cdoTi  *r*ctuar  en  «1  (U-gc»***  lento  d# 

C El  «ut«uti  >imji  un#  Ubii  *j| tiple  quD  ,n 

C ■ yucf*  [ir  un  (Ctinncr  » #rte.*rq*  rtr  drtorfif  j,C4r 

c uní  tln«*  d«  *ntr»Q#d*  * 8l, 

COMMOH  / BUFFf.fi/  ft,BtC 
CfWLEJ!  Ai4Q9ih,Bi4*j9*(,C<«iO?aj 
COMALE*  Afl<6*,4H«0  ,BBÍ  ¿4  .¿*>,1X1.44, 6*1 
EQUJ VflLElfCE  CA,Aflj  , (E^fcb»,  rc«CCf 

CHARñCTEÁ*aaa  pal 

CHAfiAC  TEH  *fi¡>  PALABRAS  1 1 0 t 

CHAfiACTÉR*  1 3 NOMBRES (30)  „ T»  EN , NOMBRE 

'^5NCE  £ rífraPC  a - COMANDO  I * t TOCEN  , MANDATO  J 

EDUI VALCNCE  í L I N ,L  INEA  ► * i PALABRAS  , PflJL  i 

COMttfflY  r'  CDHMAND  /L I N T MANDATO 
CtWON  /cntHNDI  / LOMO I TuD,  INDICE 
IRTE SEA  DIRECCION!  30  » 

c 1?ír*ceií"r!l  **  aBl**  ■ da  i caiumdp 

3W1E0ER  REQRESOr  3*>) 

C dírtteidn  d#  r *tt' r na,  un#  v#í  quv  u ejecute 

Í !JÍ™  lO«Ittl,0I V rs  l«  . tt I RECTA  , 1 N VERSA 

INTEfiER  EfiCT tfi IR pLEEft,OPEftftCIONr SALIR, CARACTER 
tNTEDER  fi£flLOinAB1eiWJU,RECCieClL,PDLOReCpNEO 

CMlfcOER  lOSAe.LOGAia.FiEñLO.IhAlJo 

INTE6ER  CERO  , CAMBIO .CutUffl^DC.OROQ,  OFFT*  INVfFT 

c D#t#«  pin  t#«  epmt^ntt»  qus  » ocuparan  tn  *1  prc^nu, 

UATA  SUMA,  RECTA  .MULTIPLICACION -DIVISION/ i ,í  3 */ 

CATA  ESCRtfilR.LEEn/l.e/ 

DATA  SALIR/ 1 / 

PATA  REALOIrtAG.coWJUjftECOPaL+paLQREC, 
«KLDSAe^lJIMiOpW^IÁw/^b.f.t.VplÓ,  11,13,13/ 

DATA  CERO *£AM6lO,CLJADfiÁIKJ, GAGO , DFFT,  INUPFT7  i*. 15*  1*,  IT,  IB,  19/ 


C MAmíímí  íUiilWrD  d*  et?**TKÍVl  Qué  I i*n#  1*  tihii 
DATA  rtfil  )MD  /jo/ 

C Ci.jr.jrKlL.rj  pur  ion  #e»pt#dB«. 

PATA  NCJMBREfi/'&iSCO 

* 'Simo 

* 1 RESTO 
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c 

c 


‘ RESTAR 
-DIVIDIR 
ÍJtvtflG 
‘MULTIPLICO 
r!SVt  i aligar 
* NESA» 
■SALIRSE 


■RESTE 

■DIVIDE 

■NULTrecTOue 

■MOT 

■NIEGUE 

‘TERMINAR 

■ACADAft 

*EXIT 

£ND 


-CAMBIE 

■TRASLADÉ 

'TRASLADO 

■CAHblAfr 

'LtHADlTHD 


-CLEAft  % 

■SORT  ■/ 


P*l*bi-*»  qrU*  U Utillílin  #n  un*  Unii  dt  CdlHfdCf. 

DATA  PALABRAS/ 

•‘PASE  10  ftEAL-IMAO  D-] 

» ‘REAL  I MW1WR JO  CON JUQ  ■ , 

* ■ AÍKl  AÉC—POL  RECTANGULA'  , 

*'R  ROLAR  fiPPfta  OfilBEW  , 
f ORO  DER-lZO  IZQUIERDA'» 

*•  DERECHA  TRASLADE  1 , 

*#  NEGAR  PARTE  pLT  , 

**PFER  LEER  ESCRIBIR  * f 


C D 1 r *e  c 1 unil  d*l  rtgrt«u. 

DATA  REGRESO/C*  . £4 , £A  , S*  * S4 ■ B4 , £4 , 

■34,24,34 ,BA  rtAFeA  + 0A  + aAf S4  »B4  ,EH , 

. ¡Z*t  i * EA , EA  , EA  , É*  , * SA  , a*t , EA  » EA  / 

C & i r »c  c lOp-n*»  dprd*  V9y  dnpii«  d*  dtííHllflEir  *t  lultiMti. 

DATA  DIRECCION/ 

*l»e»e,Er£r£'fl»£iL3tI»3i 

»3, 4, 4,3,13,13,1.1,1^1* 

*Í,V,9.9.9,IA,lÉi,lB,eíl/ 


C COMIENZO  EL  PROCESO  CON  EL  AUTOMATA, 
30  CONTINUE 

CAll  OLRHCR 

c ********* PRUEBA  *** 

CO  ID  K*i,3 


7 l 


CALL  PATRON ih, A > 

00 

DO  I*  J-l .6* 

PRJWT  -.REALIAi  ¡aéA-JH 
14  CONTINUÉ 

FRINT  - 
1 3 C£hT  ] NltE 
13  CONTINUE 

PRINT  *,,F|M' 

STOP 

CALL  LEERLlrc  *""* 

índice* i 

CALL  6CANMCA 

juhp^i 

j-i 

ao  CONTINUE 

LUQAR-I 

fc“  ICHARl  COMANDO  I | , JH 
L-ICMAflíMWjWITOÍ  LUGAR)  1 
IF  { K.EQ.L)  THEW 

C V*r-a*Bift  «1  »tg  ■«,  ai  cea-ande  qua  t*  tnwoco 

10  CONTINUE 

tfm I CHAR I HftNDATO  í I ) t 
L- 1 CMflfi  £ C0NANDO  f I , J > j- 
IF  ÍK.CQ-LJ  LOGftft- L UOflR  +■  1 

l*t  + l 

IF  ILuGAA.EO. LONGITUD]  JUHPaDtftECCIOWt J J 
IF  ILiAJAft.EG.LQNBnuO)  00  TO  40 
IF  U.Lf  ,LQN0ITUD-H  > 60  TO  10 
ENOJF 

J*J~I 

IF  t J.LE.KÍLlíND)  OO  TO  ÉO 
m fü  lOO 
40  CONTINUE 

C 

C -# *a*l**a****+,*»**. 

c C&ntamp  at  pPQCWMkianeo  dcr]  CfcMndo 

IF  1 1 JUWP.LT.l ) .AND, (JUNF.GT.aAJ 160  Tfl  100 

SO  TQ  ( 1 00. ®>0. 300,400,300**00,-700*800,^00, 1000,  II OO. 

• IfiOOi,  1300*  1400, 1500,1*00,  i 700,  |gOO,  I OOO . SOOO , 5 1O0 * 
*5200,8300,30  !tJuF  ' * 

C nunca  perdra  rf^oiar  da  I.  «al  ida  ,1*1  n,tqa«tl.  * 

C 

*-*.*•* a***,*****#**,/ 

loo  CONTINUE 

OPERACION-^ 

CALL  FlLESlDFERACÍON.MQNBfltr.A] 

PfUrn  *,LIN 
PRINT  * 

SO  TO  30 

C 

MO  CONTINUE 

300  CONTINUE 
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400  CON! INUE 
300  CONTINUE 

c Le-»  cpMrwto»  invfreadtft  «en  20p , 300 , 400  y 500  par  a *u*a, 

<J  rtst*,niu 1 1 Ipl itaciOn  y diyiviGn, 

CALL  SCANNER 

CfittflCf  EA-  ECHAR ¡AWOftTO  f I I > 
n FP  INI  * , MANCA  TO  I 1 I ,T0*'1¡:N,  INDICE 

ü PAUSÉ  "ESTOV  EN  200  AL  500' 

If  ( CARACTER. ED. *3)  tHEN 
C Era  vi  6uff*r  A. 

CALL  SCANNEF 

U PRJNT  *,MANDATOl I 1 . INDICE .TOKEN . ichaf  Lmansl ■ te  C l i ) 

ü PAUSE  ' INGRESO  AL  TIPO  DE  OPERACION' 

□PC RACION-SAL  IR 

CfiNÁí:  í£fi>  ICHARl  MANDATOÍ  | 1 > 

IF  lCAAACTPR.EO.43l  OPERAC I ON-SUHA 
IF  i,CAftAC1€R.EQ.43)  OPEA  ACION-RESTA 
IF  LCAAACTEFt.eO.42i  OPERAC  IQN-rtuLT  IPL  ICACIOW 
IF  1 CARACTER ,EP,5Sr  OPFRACTDN-Í I VISION 
IF  i CARACTER, £0,47)  DFEftAC  I I V 1 S 1 ON 
IF  t OPERAC  E CU , EO , SAL IR í DO  TO  |O0 
C Ei* a un  aigno  d«  ú^rieien  vél  utí  ■ 

n PAfNT  *, I NI DI CE 

CALL  BCANNER 

CARACTER- I CMAR l NAW0ATÜ < 1 1 ) 

D RRINT  *,!«N0AT0lll,lWOÍCEtTacEN(0P£ltACIWI 

O PAUSE  ’DECDDÍFIOUE  LA  OPERACION  uE RENOS  DONDE  VA«OS' 

IF  E CARAC  TER , £0  , *3  ) CALL  FUNCIONE OPERACION, A ,A ) 

IF  1 CARACTER  .EO-  RAÍ  CALL  FUNd ÍQNl OPERACION, A,  tf  ) 

IF  (CARACTER. EÜ.41)  CALL  f UNC I OWt OPERAC 1 0« , A , C S 
I F i 1 CARAC  TER ,ME  + *3  I , OH  . i CARPÍ IER . NE . A* I -OR . 

* ( CARACTER. N£. 47 1 ) QQ  TO  lOO 

JUMf>*R£GR£  50 1 JUNP  > 

O PRINT  *,JUMF 

D PAUSE  "VEREtiOS  OOND€  SALTO' 

00  TQ  40 

£NDIF 

|F  ( CARACTER, EO. 44)  THEN 
C Era  mi  Burfar  B. 

CALL  5CANWER 
OPERAC ION- SAL  IR 

CABACTER-ICHÁR  ÍrtANDATOÍ  I ) 1 

IF  LCARACI6R.E0.43)  OPERAC  tCPI-SUMA 

(F  ( CARACt  EA . EO , 45  i OPERAC I QN- RESTA 

IF  1 CARACTER.  EO.42.1  QP£ftAC  IQN-tflJL  T I FL I CAC I QN 

|F  ( CARACTER. 20.58)  OPERAC I OM-O IV I 6 I ON 

IF  l CARACTER. EO. 471  OPERAC ÍON-D I V I 9 I ÜN 

|F  < OPERAC I0N.EQ.SA4.ER)  SO  TO  100 

C Era  un  alann  d«  opar-acion  vAliao, 

CALL  BCANNER 

CARACTER- ICHARl WANBATQí 1 1 1 

IF  I CARACTER. EO. *5)  CALL  FUUC ION i OPERAC t QN , D , A > 

IF  ( CARACTER . ÉQ , **  ) CALL  FUNC  I ON ( OPERAC I Oh . E) ( B) 

IF  i CARACTER.  Eli,  *7)  CÁU  FUNC  I ON  t OPERAC  t ON , B . C 1 
IF  ( ( CARACTER. NE, *3 ),QA.Í CARACTER. NE.Afc)  .GR, 
p 1 CARACTER1 . !C  , *7 1 1 QQ  TQ  IOO 

JutlP-  REGRESO  f JUMp ) 

SO  TO  40 
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EHJHF 

IF  (£ARACTER,£C.fi7»  THEN 
Er 4 «1  buffm-  6. 


Cflíí  SCAWNEft 

DFÍ PAC ION-SAL  IR 

CAR ÍC  f ER-  I CHflft  [ KANDAT-Q  f | } j, 

If  ( CARACTER. El».fi3J  OPERACION-SUMA 

IF  (CflRftCTER.EQ.451  QRERftC  1 Wt-Re 5T A 

íf  i C.ARAC  t Eft  , CQ . fi£  j OPERACION» MULTIPLICACION 

IF  ¡. CARAC TEft . EO . 50  > QPERftC ION-DI VIS Itw 

tf  ( CARACTER „EQ. 47»  OPERACION- Q¡ VISION 

*F  iOPERrtCÍON.ED.'SALlftl  GQ  TQ  lo* 

C Era  u iv  b l^Tio  d*  cp*r « im  válida. 

CALL  SCANNgR 

MANDAT0 1 1 ) J 

ÍF  (CflUACíeií.CO.feiSl  CALL  F 1JNC I Qftf t QFERAC J qn , C . A ) 
JF  * t ARAC  T E Ft ,.  £Q ..  fifi  > CALL  FUÑE  ION  (QPEftAC  3 DM  C En 
ÍF  iCARACTER.ED.fi71  CALL  FUHC ION tOPÉRÁC I DM ,C .£  1 
if  S CCAftñCTER.NE.i5>  .DR.  1 CAftACTEH.NE.ifil.Ml. 

* ( CARACTER .NÉ.fi?» 1 SO  TO  100 


JUMR-REORÉSDiJUnfrJ 

so  td  ac 

END1F 
BQ  TO  * t>0 


C 

too  LÜnF ENUE 

C Cnnanda  pár»  l»ir  *n  ttilffar*. 

CALL  BCAWNER 

□PEAAC  lOWHJEEft 

if  ÍLQNGITUP.GT.Í  í SO  TO  10* 

IF  t CARACTER.  EG.ÜI  THEN 

W»M- ■ BUFFEftrt . DAT - 

CALL  F I L£9  7 OPERfiC  1 0*1 , NOMBRE  t A ) 

JUtlP-ftCGRESOÍ  JUNA* 

OO  TO  LO 
END  JF 

IF  (CARACTER  .£4).  ¿6»  THEhj 

N0K6RE-  ' fiuFF£f,  f, . dax  • 

CACL  F ILE  E ( OftERAC  E ÜM . NOMBRE , 0 1 
JUMA-REHRESO  í JUrtP  Y 
so  to  «o 
ENDÍF 

IF  <CARACTER.EQ.i7í  THE H 

NqHBftE  - - BUFFEAJT  . DAT1  r 

CALL  F IL  ES  1 ÜRERAC I OH , WOflBftE  . C 7 

JUHF-REGRESÜ 1 JUHP  j 

60  TO  fi* 

ÉNDIF 

C Error  fin  Ucturi  d«  burrera, 

&0  TD  100 


® Eomtuía  P*r*  HCrlblr  le*  buTrcra  *n  di  acá. 

TOO  CONTINUE 

CALL  SCflífiiER 
OPEPAC I ON-ESCR I B I A 
IF  (LONajTUD.BT.  ! J (50  TO  1*0 
IF  ICñftACiÉR.Eq.fiSf  THEN 

HOMBRE»  k Buf  Fffi/L . MT  * 
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£ 

£ 

£ 


C«_1_  FlLE9IClRERACI0H.N0WEiAt.AJ 
J UMPHR&BUOD  ( JUNP  i 
t.Ü  TÚ  40 
EN&1F 

IF  [CARAC7EFt.E0.ifi»  r HLK 

NOrtBRE- ■ BUFF£R6 . DAT ■ 

CALL  F J LES  í OPERAD ION , NOMBRE ,B) 
JL^P-REOftE  BD 1 JUMA 1 
BO  TO  4* 

EHOIF 

If  (CARACTER. EG. 67)  T*£H 

NOMBRE» - BUFFERC . OAt  * 

CALL  F I UE& i QPERAC I Otí , lAínSRE  ,C  1 
JUNP-REQRESO ( JUMR > 

OO  TO  40 
ENOIF 


A cartiruáUdfi  tlttn  %tr dtl  le*  fiintieflM  HU*  »»  p&ilblr 
h*b*r  coa  loa  bMff-r*,  BUHA,  RESTA,  Diyisiew  y NUL T 1 FL ICAC ItW . 
CONTINUE 
CALL  5CANNER 
IF  [ CARACTER .E O - *3  » THEN 
CALL 

□reraciow-baltr 

IF  < CARACTER. EO. 43»  QRER0C I OH- 3UHA 

ÍF  íCAAA£TiR,ea.*5>  OPEPAC I ON-REÍÍTA 

IF  < CARACTER. EO.fiCl  OFÉftAC JON*HULT I FL t CACI ON 

IF  i. CARAC TER . EO  + 50  » DPÉRAC  EOft-03  * 151 QN 

IF  1 CARACTER. EO. 471  M*CRAC I DN-D I Vi 61 Oh 

IF  (OPERAClÜN.EO.añLIR'i  60  TÜ  1*0 

CALL 

|F  ( CARACTER  - E'Q  ■ ¿5- í CALL  FUKUORlSKMC ICW.  A.  A) 

IF  [CARACTER, ES. fiA>  CALL  F U«C I ON ( OPÉÁÁC ION. A. 0> 

!F  ( CARACTER .EB-Í7Í  CALL  FUNtlQKtQFÉRAC IQN,ATC1 
\f  t (CARACTER .NE .fi5 3 .OR. i CARACTER . NE . ¿¿ » .OR . 

4 (GflRACTEft.NE.i7t  i QO  TQ  lOÍ» 

JUHR-REeRteOiJUHRt 


ENOIF 

IF  ( CARACTER. Efl.iit  THEN 
EfltL  SCANNER 
OPEPAC (OH-BALIP 

IF  1 CARACTER. EQ. 43 J QFERAC I «OEÜrtfl 
|F  .CAftACTER.E0.45»  QPERAC IW-flfB í A 
|F  i CARACTER  . EO . 4-3 1 QPEftflC  lON*(1ULT  I RL I EAG I ON 
IF  (CARACTER.Ea.5Bi  QF^fiRCÍOtT-OlVlSlON 
IF  (CARACTER. £0,47»  OPEftAC 1 DN-D I VISION. 

IF  (0PERrtCiaN.EQ.3AL  IR J SO  ÍO  1*0 
^ ^CAf#4ER 

IF  í CARACTER. EO.fiSl  CALL  FU4C 7 DN í OFERAC I ON , B, Al 
IF  tCflftACTER.EQ.ii»  CALL  FUNCIQNtOlPfiRftC  ION.  B.B  1 
|F  [ C AAfiC 1 ER . EO . i7 7 CALL  FgNCI0WtnPEPAClQN.6,C i 
If  1 ( CARACTER. NE-fiS» .OR. [ CARACTER. ME. 66 } .OR. 

» [CfiPAf:TCP,í(E,i7H  00  TO  1*0 

JunP^PEQREaO ( JUNP ) 
ti*  TQ  4* 


ENDIF 

Ü ( CftfiAC Tift .EO , t.7  I THEN 


QPERrtC ION-BAL  IR 

IF  (CAftACTÉR*É0*43  » QPERAO ON^SUHA 
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IF  ÉCAftflCTEFr.EÍKi^J.  OPERACION-RESTA 
JF  (CARACTER.  EQ**t£>  ÜPE  PAC  I Un^mul  t 1 Pí_  ! C AC  I Üti 
IF  l CARACTER,  EO,  S0 ) OPERAC  I ON"  D 1 L I S I DN 
IF  (CARACTEft.£CL*7J  OPERAC ION-* (VISION 
IF  (OPERACION rE0. SALI Rt  GO  TO  109 
CALL  5CANNCR 

IF  (CARACTER ..£0.65  I CALL  FUÑÉ ION  (OPERAC lOK^t, Al 
IF  I CARAOT £R  r EO . 6* > CñLL  FUNCION  I C)P£R*C  ION  ,£.  B > 
IF  í CARACTER. EB. 47 > CALL  FÜNC I E3KÍ  DPERAC í ON , C . C ) 
IF  < < CARACTER.  NE  + 65  KOft.fCARACTER.fC,  *61  .OR. 

(CASACTER.NE.471 1 00  TO  HXi 
J ulF"  ft£  ORE  SO  ( JUMP  I 
BC  TO  40 


Precedí ai ante  para  obtener  filtra  »n  imagen  p bien  FFT  directa 


C 

900 


LA.  Adame  cambie*  vario*  en  úna  «atril* 


Intercambie  parte  I 

DFERAC I W-AEAL01  MAC 

CALL  BOF“ 

GO  TO  l 


. par  IrtASÍNARIA 


C Pr peed i mi ente  para  obtener  *1  conjugado  da  una  aatris* 

lOOO  CONTINUE 

OPERAC ÍDW-CQNJU 

SO  TO  9 

c 

C Rreced  luden  te  para  obtener  a]  cambio  da  rectangular  a polar* 
1 100  QOMTENUEf 

□FERAC  lOff-RECOPOL 
go  rq  i 
C 

C Procadi,» tente  para  convartir  da  pelar  a rae tangid ar  . 

IfiO*  CONTINUE 

OPERAC  IQN-PGLOREC 

□O  TO  | 

C 

£ Precedía!  lento  para  ™g»r  urna  matriz* 

1300  CONTINUE 

ikwcioi  tma 

QO  TÜ  1 


C Procedimiento  para  «acar  al  logaritmo  en  fai»  E dr  una  matrla 

C Procedimiento  para  eb tenar  ai  logaritmo  ar  pasa  lo  d*  una  asir 

It. 

14*0  CONTINUE 
1500  CONTINUE 

CALL  SCAMNER 
► - INDEX  ? PAL  , TCIlEN  S 
IF  IK.EQ.Il  smLn 

OPERAC 1 Qft-SAL  IR 
CALL  SCAMNER 
P>  INDEX  1 PAL  ,TQFEN1 
IF  fF.Eú.4 i OPERAC ION- LOGRE 
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c 

c 

c 

140* 

170* 


IF  ll:  + tÜ**J  OPERAC  ION-LOBA  10 
I F I OPERAC  ION *EE * BAL I R > GQ  TO  lGO 
m TO  1 

EL  SE  ÍF  (*:*£&*  4 t THCN 

OPERACION-LOGRE 
DO  TQ  t 
EL3E1FCF.E0.61  THEN 

OPERAC  [Ot-LfflSAlú 
QO  TD  1 

ENDIF 
QO  TQ  lOO 


Procedí  miarte  para  ¡lavar  a tare  la  parte  REAL  d ¡>  una  matr 

b La  pert*  (HAG 

CONTINUE 

CONTINUE 

CALL  SCANNER 

K-XNSCXLFttLtTOKeO 

IF  CP ,ED, liu i THEN 

CALL  SCANNER 
□PERAC ION-SALIR 
K- I NOE  X t PAL . TOkEN  í 

íf  ÍK , EQ* I 1 1 OPERACION-REAlD 
IPCF.EO.il>  DPERAC I OH— I MAGO 
IFf SALIR -EO.OPERftCIDmSa  TO  IOQ 
GO  TQ  1 
EL5E t F í K * E* ,111  THEN 

OPERAC  ION- R£A(-0 
GO  TO  I 
ELS£  I F r F EO * 1 A ) THEN 

OPERAC ION- JUAGO 
GO  TÚ  1 

EMOIf 


Precedí  «i  *nta  para  barrer  el  buffer  indi cade. 
CONT INUE 
CALL  SCAMNER 

K-lNOCitiPALi  toben  í 

lFíH*E&.5Si  r«N 

OPERAC  ¡CHM-CCRÚ 
GO  TO  l 

EL SE 1 F t K , EO+ 1 53 ) THEN 
QO  TO  l 
ELSEIFÍK.ÉÚ. I S TREN 

INDI  CE*  INDICE-I 
GO  Id  I 
E L SE  1 F C > . . CO  * 2 l t&H 

IND1CE-ÍNPICF-1 
GO  TO  | 

EtSEIFCK.EQ.TM  TH£N 

IwOlcC"  iNtil  CE-1 
SO  TQ  1 

ENDIF 
SO  TO  tOO 


e Procedimiento  pata  Intercambiar  parte  derecba  per 

£ la  izquierda  en  una  imagen. 

|W0  CONTINUE 

ÜF'ERAC  I DW-CANB1  O 
QD  TÜ  I 
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zoos 


Precedí  bu  wntP-  p*r*  **e.«r  *¡  cu«ír*£ífr  dv  L*  m«DTUtud. 

CONTINUE 

□PERAC  I Om-CUABRADO 
0.0  TO  1 


C Pr  Dcsduutnti.  pif*  pe>mr  *n  eert>  mí  cuaen  de  Ié  Iitigcn. 

eiOÍ  CONTINUE 

OPEJWCI0N-0HBO 
QC  T0  | 


C Pi'ccadtaitntD  per*  &bt*n«r  i t trerVAforAed* 

C ém  FOURIER  d*  1* 
eaao  continué 

0PERACIQN-DFF1 
GO  TO  1 


C Procedí*] *nt o pare  obtener  I*  trinferiudi  Invtfiá  de  füurjen. 

£300  CONTINUE 

OPERACION»  INVFFT 
GO  TO  I 


Prficiauimto  ptr*  procíMf  todo*  lee  cDnjndci  indicidet. 
CONTINUE 

sríJKSU,,, « ™ ^ 

IF  t CARACTER  .EQ.áOii  CAI_L  FUNC I OM 1DPERAC I Sn,a,B) 

IF  í CARACTER.  ED.  ói)  CALI  FLINC  I QRt  ÜPERAC  I TJN , E>  + fl  > 

IF  {CARACTERES,  b?)  CALL  FlMCIONtOPERACION,!:.  Al 
IF  i ¡CARACTER.NÉ.éS)  .Oft.i  CARACTER  .n£  .6*  ) .«*. 

* (CARACTER.  #€.*?**  THEN 

QC  TQ  100 
ELBE 

juMc-ftEaftEsn  t juw  i 

GO  TO  -G 


Jump-i 

GO  TO  'rO 

EMD 

€ND 


GUSROUT1NE  FUnC  1 QNí  COD I 5D  - A t B > 

C 

c 

C ESTA  EUBRÜUT  JnE  EFECTUA  LAB  DPERACJNES  AR  f TM6T I CAS 

C SOBRE  EL  PLANO  tliABEN  DE  MATRICES  A V p ANSAS  COMPLEJAS, 

C CdDTOO*  l 

C SE  EFECTUA  LA  SLH1*  DE  NATA  ICES  R-A  * S 

C 

c coQieo-e 

C SE  EFECTUA  LA  RESTA  0E  MATRICES  A-A  - S 
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A-A  * Ej 


C COOlBO-3 

C SE  EFECTUA  LA  MULTIPLICACION  DE  MATRICES 

C 

C COTIOO-* 

C SE  EFECTUA  Lfi  DIVISION  DE.  MATRICES  A-A  í B 

C 

c CODIGO- 3 

C S£  EFECTUA  EL  TRASLADO  DE  A — > P y S “>  A 

C 

C CtttHBQ-4 

C SE  EFECTUA  EL  TRASLADO  DE  PARTE  REAL  POR  IMAGINARIA 

C 

C CODIÓO-7 

C BE  EFECTUA  EL  CÁK&IÚ  DE  MATRIZ  POR  EL  CONJUGADO, 

C 

C CODIGO- B 

C SE  EFECTUA  EL  CAHBTO  DE  RECTANGULAR  A POLAR. 

C 

C CODIGO- 9 

C SE  EFECTUA  EL  CAMBIO  DE  F13LAR  A RECTANGULAR, 

C 

C CODIGO- i<s 

C SE  EFECTUA  EL  NEGADO  DE  LA  MATRIZ. 

C 

c CODIGO-I  1 

C SE  EFECTUA  EL  LOGARITMO  EN  BASE  I. 

C 

C COOIGQ-ia 

C SE  EFECTUA  EL  LOGARITMO  EN  BASE  i O. 

C 

C CODIGO- 13 

C SE  EFECTUA  LA  PUESTA  ñ CERO  DE  LA  PARTE  REAL. 

c 

C CODIGO- l* 

C S£  EFECTUA  LA  PUESTA  A CERO  DE  LA  PARTE  IMAGINARIA 

C 

C CODIGO- IS 

C SE  BORRA  EL.  BUFFER. 

C 

C COBUÜO-li 

C SE  INTERCAMBIA  PARTS  DERECHA  POR  IZQUIERDA  EN  LA  IMAGEN. 

C 

c cooíeo-iT 

C SE  SACA  EL  CUftDfiApíl  {£  LA  MAGNITUD. 

C 

C coDitiO-ia 

c se  PONE  A CERO  EL  ÜMGEN. 

C 

CÜtWtGM  /COMMAMD/L1NEA.C0MATÚ3Q 
COHPLEX  RESPÍ40VEI  .AAt^A.pi.)  .BBiM.frEl 
CGMPlE*  mWdi.BHüTti 
EOUIUALSMCÉ  <RESP+AA,tÍD 
INTE SER  CODIGO 
C 

lfiC001G0.ED.il  THEN 
C SUMA  DE  MATRICES  COMPLEJAS 

m LO  1-1  .-«m.J 

Atll-Aííl  * BU  1 

1G  CONTINUE 
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ELSCIf  ÍODDIGQ.&^SJ  then 
C fiESÍft  EJE  MATRICES 

DO  SO  1-1,4096,1 
Aíl l-Al  ü - B£I> 

20  CONTINUE 

ÉL5E1F  tCGDISa.ÉQ.31  THEN 
C MutTlPLIC ACION  D€  MATRICES 

DD  30  1-1,4076,1 
A£ Il-Ai í I * B( II 

30  CONT IMJE 

KLSÉEF  LCaDIOa.EO.43  TMEN 
C 0 IV IS  ION  DE  MATRICES 

DO  40  1-1,4096,1 
Atl  t-AÍJ  ¡l  / BU) 

90  CONTINUE 

ELEES F <CQDlBa,Ea,5>  ThEh 
C BE  HACE  EL  [ NTERCAffB  I O K A — J B Y B — : > A 

DD  50  |-í ,4074, l 
RE3PS Il-AIEI 
flU>-6í  II 
Bí I l-REBPí li 

30  CONTINUE 

BETURW 

ELSEIF  (CODIQO.EO.Bf  TMÉN 
C INTERCAMBIO  PARTE  REAL  POR  I HABI  MAR  la 

PRINT  1 nturcjatlí  o*  Hirl.  REAL  per  1IWB1NWUA  * 

OO  101  1-1,4076 
PFt£aL-REAL£A(I  I I 

pinAe-Aiiwi6iA«m 
Al  13-CMPL¡H;FIMAG,PREAL3 
IO|  CONTINUE 
RETuRn 

EL9EIB  (CDOlBCUEGUTl  THEN 
G CANSIO  LA  MATRIZ  PoP  EL  CONJUGADO, 
eoo  Pftlm  m „ • CimfaiO  Iftn;  per  «I  COMJUGADG' 

M EOl  1-1 ,407fc 
A( I )-COMJQ<ñ<  333 
SOI  CONTINUE 
RETURN 

ELSEIF  ICOOieO.EÓ.BI  THEN 
C CAMBIO  OE  ftíttANtíULAP  A POLA», 

300  PR  [NT  i , * CiMta  di  rittinguUi'  » polar  " 

BO  30 í 1-1,407* 

WAL0A-GAB51 A lili 

AHGUL0-9O.OafilGH£  Ua#«lNAG(A(l  1 I I 
RPÉAL-PEALSAS Ilí 
PlrtAG-AÍMAOlAl 13 1 

IF  [REALÍAIIJJ ,NE,0,0í  ASM0ULO-37 . ETSTS*  ATANS  i P E MAS  . FH*f  AL  » 
Al  I 3 “CMPL K 1 VALOR , ANÜULÚ I 

301  CONTINUE 
RETLMN 

ELSEIF  ICÜDICO.EO.7?  THEN 
C CAMBID  DE  POLAR  A RECTANGULAR. 

AOÓ  PtrtMI  </  Comlirta  de  p*I*r  A rectangular  1 
DO  401  I* ],4096 
AW3ULO-Q.Ú17453¿9*ñ|nAG4  AS  I > J 
PREAL*CDfi£flMGlJLO]  ettEALtflt  I 1 I 
P IttAG-SI  NS  ANÜULÜ  > -RÉA4  SAI  III 
A S 1 1 -CMPL  X ■£  PB£  Al,  , A I NAO  | 

401  CONTINUE 


ao 


RSHJftN 

ELSÉ1F  £ Cop  ISO  +ED  + 1 0 > TIC* 

C NEGADO  DE  LA  MSATRl  f . 

500  PRINT  El  nogüila  4«  1*  ««trii  * 

DO  SOI  t - I , 407  b 

ATI 1»-A(T1 

501  CONTINUE 
RETUFrM 

ELSEIF  (C0DI00.ED.1U  RÉH 
C LOGARITMO  EN  BASÉ  É 

600  PRINT  Lagar  lima  #n  bu»  E ' 

DD  fiül  1-1,4076 

VALOR-REAL S A(I  1 I 

IF  iVALOHpEQ.O.Ol  VALOR- 1 0 . OE SO 
VALOR-ALOBS  VALOR  I 
A í 1 i -CMPL  1 1 VAL43R  ,Ó . 9 1 

601  CONTINUE 
RETUR» 

EL SE  |F  S CODIGO. EO. 1E ) THEN 

e logaritmo  en  base  10. 

7W  PÑINÍ  Logaritmo  en  til*  IO- 

DO 701  1-1,4096 
VALOft-R£AL t Aí 111 

IF  4 VALOB. £0,0 ,01  UALOR-IO. OE-PO 
VALfft-E  . 30S39  -ALOE  í VALOR  I 
Al  1 I -CMPL J< VALOR *0, Oí 

701  COMI  3 hué 

RETUNN 

ELSfc  I F S COO-I  SO  i ED,  t "3  > THEN 
C RONCA  a CERO  LA  PARTE  REAL. 

300  PRINT  poner  a c*rí  la  parta  r*al ' 

DO  SOI  1*1, 4076 
VALOF-AIHAQ1A1 III 
Allí -CMPL 1(0.0, VALOR  3 
SOI  CONTINUE 

RÉTLIfiN 

ELBEIF  ÍCDDIH0.ED.14S  THEH 
C PONER  A CEPO  LA  PARTE  IMAGINARIA. 

TOO  PRINT  •/  Parar  a caro  la  parta  ttaaqlnar  i*T 
00  Rol  1-1,4076 
VALOR-A I MAG t A T I >1 
Al  I 3 -CMELJt  < VALOR  ,0.01 
701  CONTINUE 

RETURH 

ELSEIF  SCUDIGO.EÜ.15I  PHEN 
C BORRAR  EL  BUFFLR. 

1000  PftlNT  *♦'  Barra  al  buffar  ■ 

DD  IOOI  1-1.4076 

Al.  1 .1  - CMPL3;  1 O .0,0 .0  I 

1001  GOHTINUE 
RETLffíN 

ELSEIF  fGQDlBO.Ea. 161  THEN 

C INTEPCAffBIO  PARTE  DEKRCHA  PÜR  IZQUIERDA  OE  LA  (haOCH. 

ÉOOO  FFíínT  *,j  ln ver cambia  parta  daracAa'* 

• ‘ par  par  ta  uqniP  d*  an  la  miagan r 
DD  eOOE  1-1,4096 

R£gPl  n-*s  n 

eooe  CONTINUE 

00  £<M.n  1-1,30 


(1.63-1 

no  soiH 
BAk-AAiJ.I  I 
AAU  , II-AA<JrltI 
rtñií^t  n-Bfln: 

?«i  continué 

t) O £i>P3  I»1,,40VA 
Pii  I I-RESPÍ1Í 
£003  CONTINUE: 

FtETUfiN 

ELSEIf  iCOEJIGü.EO,  171  THEW 
C CUADRADO  DE  LA  HAÜNI  TlhÜ  . 

30t>0  PfftlifT  *,’  CiiÉcfr (Míe  -da  ] j magnitud  * 

00  3001  I*1,LQVA 
mí  I-CAES  (Ai  1 > )**e 
3001  CONTINUE 
RETURN 

ELSEIF  <CnP[GÜ.EÜ. L0>  ÍHtN 
C PONER  A CERO  EL  ORIGEN  DÉ  LA  [MACEN, 

4000  PRÍNT  Punir  a Cüro  fll  Gflqcn  " 

ACLORYI.iO.O.O.OI 

RETUHN 

ENDIF 

ENQ 

BUBROUTINC  PRTftOmPRI  „A) 

e esta  rutina  genera  en  la  matriz  -a-  un  patrón  be  test  i 

C 

C PAT  - Í — ? FRANJAS  HORIZONTALES 

C 

C PAT  - E **?  FRANJAS  VERTICaJES 

C 

C PAT  - 3 a*  AJEDREZ 

1RTESEA  PAT, COMÍA 
REAL  RL  ANCO 

REAL  negro 
REAL  VALOR, I«A 
COMPLEX  A(AOÍO> 

blanco  - IO*S 

NEBRO  - -10,0 

if  <PAT,eo,ij  tmen 

C PATRON  Ü£  FRANJAS  HORIZONTALES 

CONTADO 
valor- Blanco 
00  10  l-lT4P0é,l 
COKTA-CÜHTA* I 
1F  (CÜNTA .EO ,256  r TMEN 
CONTA-Q 
VALOR  - “VALOR 
ENDIF 

IMA-AItlñQj  Al  IT* 

Al  I l-CMFLX (VALOR, 1 NA > 

10  CONTINUE 

ELSEIF  '.PAT.ÉO.ej  1MEN 
C PATRON  DÉ  FRANJAS  VERTICALES 

CONTA-O 
VALOR-BLANCO 


S2 
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00  £0  tM**ÜS<t.Í 
COHÍA-tONTA+l 
IF  1 CÜNTA  lEO  , A ? THEn 
COMTA-Q 
VALOR  - -VALOR 
ENDIF 

IMA-AIMAGIAIISI 
A(  i)  -CHfa,XíV«,OR:tIMAl 
SO  CONTINUE 

ELSEIF  I PAT  . E O . 3 > Tm£N 
C PATRON-  ü€  A.JE0R£2 

DO  30  l-ltAOflt,l 

1P  I < I HOB  £f,EQ,U  THEN 

IF  ÍÍF1LA  NOD  £S ,E0. I I THEN 
VALOR  - blanco 

ElSE 

valor  - m£0RÚ 

ENDIF 

ELSEIF  HFILA  MOD  Si  ,É0.  I J ÍHÉN 
VALOR  - NEGRO 

ELSE 

VALOR  - BLANCO 

ÉNOIF 

A( I I -CMPLF  C VALOR, VALOR] 

30  CONTINUE 

ene  i f 

RETURN 

ÉNfi 

C 


SUBRÜUTINE  F I LEE ( OPGRAC I ON , ARCHIVO , MATRIZ > 

ESTA  RUTINA  HACE  LA  APERTURA,  CERRADO*  LECTURA 
□ ESCRITURA  de  UN  ARCHIVO.  LA  LECTURA  □ ESCRITURA 
SE  HACE  POR  MEDIO  PE  LA  MATRIZ  B, 


OPERACION- 1 

ENTONCES  SE  ABRE  UN  ARCHIVO. 

OPERACIOW-e 

ÉNTCHCCS  SE  CIERRA  UN  ARCHIVO. 

PPCPAClÜN-3 

ENTONTES  SÉ  LEE  UN  ARCHIVO. 


OPEMCialM 

ENTONCES  Sfc  ESCRIBE  UN  ARCHIVO, 

OPERACIDN-3 

ENTONCES  SE  ESCRIBE  UN  MENSAJE 
DÉ  COMANDO  MAL  F ORNADO . 

LOS  DATOS  LEIDOS  Q HA  ESCRIBIR  BE  DEBEN 
ENTREGAR  TTN  LA  MATRIZ  B DEL 

BuFFER  arch. 

• r « «••*!•* av ■** ** .*»###** *»»*#»»•»»»»*** ** 

COMPLEÍ  B I 409&1 * BB  t *A , ¿h 1 * MATR I Z i AOVA  J 
EQUIVALEN^  i f,  n (¡6  J 


CHARflCTER*l£  ARCHIVO 
INTEGER  OPERACION 


EF  ¡OPERACEOW.EQ.I 5 THÉN 
ABRO  LHK  flflCMWD, 

□PEN i UN  n -4 ile-arch I va . 

• ACCESS- ■ DI RECT 1 , AECL-3 12 > 

RETORN 

EL  SE  ÍP  I DF-EAftC  J OW . EÚ  } THEN 

CERRAR  UN  ARCHIVO. 

CLQSE<UNIT-«it 

RETÜRN 

EU5EIF  <0PERACI0N.ED.3!  f«N 
> ARCHIVO  EN  LA  IWRI2  DEL  SUFFEft  AíiCH  LLANADA  6, 
DO  ID  1*1)64,1 
UO  10  J-l.64,1 
ftEADf UNI T-A , REC- I )EB( I, Jl 
CONTINUE 
DO  1 I I-I  ,40-?* 

MflTRIZí  | í-BE  | | 

CONTINUE 


SE  LEE  I 


ELSEIF  (OP£RACION.EQ.«h  Thcn¡ 
ESCRIBO  EL  &UFFER  AFkX  EN  EL  DISCO . 

DO  El  I-1.4OT6 

bu  i-nAmm  n 

CONTE  hLIE 
DO  EO 

&0  £0  J-l)64,l 

l*IR  I TE  C UN  j 1*4,  REC*  I J bp  I I ^ J I 

CONTINUE 

WETUftN 

_ EL  BE I F 4DFERACIQM.E0.51  T«EN 

EBCRI60  UNA  PANTALLA  DICIENDO  QUE  EL  EflfiQR  SE 
DEBC  A ILEGAL  COMANDO. 

CAU.  CLRSCR 
PRINT  • 

qu»  invoco  11  ltEgj|,r 

RCTURN 
ENDÍF 

RtTURN 
ÉNO 


BUBRÜUT In£  CLRSCR 
CHARACTEAvl  FILA  LIO 
CHARACTEft^lO  BORRAR 
EOUJVALENCE  < F I LA . BORRAR  i 

DATA  FU.A/í-SC't7-  1 B%Z-5»-  ,Z'3e-  ,2 '4A-  * 
■ Z'0O*,£-QQ%Z*0C*«2-ÜO'  ,2*00*/ 

PftlNT  IGO.,  BORRAR 
l«  FDftMATífllO 3 
EMU 

C ***** 


SUBHOUT I «E  GQTOJf  Y U , J , PO«Éft  1 
CMARACTERu  FILAIS* 

CHAftñC  TER*£r  F 1LA1  .FILA? 
CHARAC  TEft*  9 PONER 
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OATA  F ILA i 1 I , FILALE  3 ,F ILAI 3) /Z *SC ' „ Z’ t © ’ ,2 '5B' / 

lt-1-  [|/iOit]0  4 ICHARÍ'0'I 

F ] LA1 4 C HAR { I y i Q + ICHAftt'O'TliVCHAfilKI 

*>J-  tJ/ 103*10  4 ICNARí'O'l 

F I lAZ - E MAR 4 J / 1 0 * I EHAA I ' O ' tí/ /CHAfl f K ) 

PGNER-FÍLAÍ 1 J / Y F I L A í 2 h / / F I LA  t 3 3 / YF  E LA I ^ Y " t ’ f f¥ILAS//‘H* 

RETURN 

ENO 


5UBRERJTINE  GOTO 
£HAftftCTER*l  FILAílll 
CHARACTEFl* I I PONER 
EOulvALENCE  l F I LA , PONER i 

DATA  F lE-A/Z’SC'  rZr  IB'  ,Z"  3B  ■ *Z  *30'  «2*36  ’ .Z*3B* ,Zm  31  - .Z'SE’  * 

* Z'LH*  *Z*M'fZ’0O*Y 
PflíNT  E'-"f  . FlSMtP 
IOO  FORHATUh+^l 
EM& 

SJSPQlJt  I NF  FftAHE  i l***eftLEFT  ff , UPPERLEF  TV  * LOWERfi  IGMTX  ,LONERR  I OH  Y í 
C 

C PARAHE  IROS  PARA  MANEJAR  EL  CUfiSOflr 
E UPPERLEFTlí, 

C EJE  IZOuIPWKl  i 

O 

C UPPERUEFTVi 

0 EJE  IZQUIEftl»  V 

c 

C LOwCRRiertKi 

C EJE  DERECHO’  PARTE  DE  ABAJO  EN  X 

C 

C L OHERfi I BHY t 

C EJE  OEJSECHO  PARTE  DE  ABAJO  EN  Y 

C 

c 

C 

ENTEGEN  UPWERLÉFT K ^UPPERLÉFT V 

CMAftACTER*!  LEFTCORNER,RTGNTCORHEL  »LINEA,fHJRALLA 
CMARACTER«9  PONER 

CftLL  OOTOiV  t UPPERLEFT )( , UPPESLEFT  Y ,Pt3NER  h 

linea-íh- 

leftcoricl-ih. 

RIOHTCfinNEL-lH. 

MURALLA- IH» 

PRINT  |00*¿£FTCORNEL 

00  10  I-UPPERLEFTX4-I  , lOwERR  I GHTJ!  - 1 1 1 
PRINÍ  IOO, LINEA 

10  CONTINUE 

TOO  FOfiHATí lH+t0OAl ü 

FfílNT  100,RIQMTCORrl£L 
m SO  I ■LPPERLEFTV+1 ^ LONERR [QMTV-1 , L 
CALL  SOTO (y l UPPERLEF  T 1 . I , POWER  f 
FRINt  |iXi,t*UftALLÁ 
CALL  GOTO  KY  i LONCRP  í 6hT  K , | , PfJNCjí  1 
PRINT  IOQ  ./MURALLA 
SO  CONTINUE 

CALL  BOTO I Y I UPPERLEF T X rLOWCRR I SHTV  t PONER > 

LEFTCDRNEL-IHL 

ftlEHTCDRNEL-lHJ. 
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PRlFíT  IOQ..LEFTC0RNCL 
M 30  I “UPPEfiLEF  T * k , LOUCRR I SKI  |f- 1 1 
«ÍINT  LOO<L|H6A 
30  CONTINUE 

PfilNT  1 00  t R t SHtCCNRtlÉL 

MIUftM 

End 


HJSRDUTINE  LÍFRl  [Fi£ 

E*t*  rulinp  lt*  un  Itnn  **  co-m*ndc  d**d* 
■ 1 (irnuul  TTY  CfijMi  UxliH  1*«  SO 
cincttrtl  fcn  L ] NEm  . 

CDMhQN  /CQHHfUf  O/i.  th 
CHflBACTCft*SO  LIN 

C 

fi£AD<*,100JLIN 
loo  pQRfiftrTíwo] 

EW 


BUBRDUTINE  SCAfliCR 
C 

C E*t*  lubfi-utin#  h«e«  mi  *******  úm  ]«  Ikn •* 

C intffedutitfj.  por  la  TTY  }u*qo  U urDítdB  A 

C *f*clu«r  *1  pmrmm  de  1*  iitM  i ir**. 

C 

COfIMON  /CQHflftND^LINEA.eaNANMS 
O0WWN  ^CWKflMOl/LQMGlTUn^tNfilCt 
CMflfMC  TER*  1 3 COHftNDÜ 

CHflBACTEfi-  í TDKEN  t 1 3 1 , L 1 MEA  1 B0  > tCR „ TABf  BLANCO 
INTEGEB  CARACTER 
IQÚIVALENCE  i TÍJKEN  * COHÁNtH)  3 
LOQICAL  fu® 

tmTft  CR.EUM'ICOyTA&.'ZOÜ'  ÉZ*«**/ 

C 

FlAtw.FHu.se. 

LOMQITUP-0 
DO  lO  1-1 , & 

10  TOfcEIHl  ti -BLANCO 

C I ijnrr'- f líi  blincoi  *1  cniílfrtío,  bl<ncm  y t»bi 

I-TWICC 

J-l 

100  CONTINUE 

CARAC  T£r-  J OMR I L INEA  í I t ) 

ÍF  HCARftCTEfi.E0.32i  .Oft,  i CARACTER . EQ,  0 i ) Tk*M 
t-I+I 
tío  to  ICK» 

ENDtF 

ROO  CONTINUE 

ÍF  (<ICHftñTUNeflUlKBE.y7)1ANDHH>lAR<LINEA!l>Í.LE1lí5>} 
THEN 

C Cattvkfertc  t$m  »|rtdKj  ] «ti  * mi yutculis  y 

€ ftonflft  dr;™s»  d*bi  nfir. 

TOFEN  [ ] ) -CHAA  i JtHftRt  LtWAC¡l  J - 3£j 

J-J+Í 
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60  TO  300 

enqíf 

c V*r  if  ic&  *t  ü*  M¥Ú*tiil*.  _ 

IF  H TCMflRfLTWtAÍ  t M .OS.fcSl +ANC- t ICMSBíLlMPA<  1 1 J»LE*^OJ  1 THEN 

C E r*  ffi^vumcu  L4  fluirdíi  *n  luo*r  **SLiro. 

T»,-ENiJi-Uf*A(  I S 
1*1*1 
>Jt| 

SO  TO  200 

END1F 

C No  ir4  I#tr«  MVúlCula  nk.  *¡nu»wH 

c ^"'"^TtCHAfiii.íwEAtl  * HBE.ASJ  .AND.  I ICHAftlLlNEftU  IHLE.StH  TUiN 

c Er«  urt  fWuwrp.  16  <Ju*rc Jo  ort  tuq#r  »qur  O. 

ICfr.ENí  Jl -LINEA II  t 
1-1+1 
J-J+l 
L-Ü  TO  E0O 
ENDIF 

c itoFfini  ou*  co*ft  *xtraft«  pu«d»  **r  ■ 

ÍF  i i t.  ICHAfi  CUMEAÍ  I i i .E0.S1  > .DB.I  ICNARtLílCAÍ  J 1 > ■OR' 

. I ICHftfitLINEftrl  J . lED.43)  .C»,  l ICHflfltL  l«A 1 1 M .ED^?t 

■..□R.  [1CHAR(Llí«flí  I J 1 .ED.í)3l!  i .»I0  . i J + EO . 1 S I tHCN 
C Es  Ljft  ligdD  6 bifn  '*  _ 

TOKEHC  Jl-LlítfAl  I S 

lHDtCS-I 
LOfífll  TUD— 3 

«eruRN 

EKJIP 

C Crt  un  aflminfltír. 

INDICE- I 
l_DNQITUt)-í-l 
RETURH 
ENO 


auBROUTIHE  FILTROBíFlUTRDÍ 

ESTA  SUPftOUTINA  CHUCULA  DI FEftEATES  TIFOS  DE 
FILTROS  P#*A  64-  rtANÉ^O  DE  IMASEn 
ÉL  UflLDfi  DC  LA  UARl A#LF  FILTRO  NOS  INDICA 
CUAL  SERA  EL  FILTRO  A PEÑERAR 

FILTRO  ,ENTDNCeg  gE  cftLCm_ft  FILTRO  SUTTtmJORTH. 

F ILTRO-e 

ENTONCES  se  calcula  EL  FILTRO  Él  ponencia*.. 

F ILTRO-a 

Entonces  se  CALCULA  EL  FILTRO  TRAPEZOIDAL . 


FILTRO** 

ENTONCES  SE  CALCULA  EL  FILTRO  LftPVAC IANO. 


COWiOW  /lUFFÉfí^  A 

CONPLEX  flf*0**> .Aac**,**} 
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EQUIUAlENCE  ( A , AEl  j 

INTE&EFcn?  PILTRO 

AOUI  DkRe  IR  UN  IF  THEN  ElJJ£  . 


CftLL  CE.RSCR 

cflLL  Qútaxr(erioTf>oNERj 

NfilTE<-,iooi 

,W  FOBHATItH*, ‘D**»  *1  vil0r  d*  Ib*  punto*  d*  certa  *rt  aja*  | ( y 

REfim  * , SOO ) COA  I C v , COR  TE  Y 
FOfmAT<£f=lü,0) 

1F  i CORTÉ í,t £,&.*>  CQRTFK-i.O 

if  ECOñTrr.LE.o.oi  come  y- 1.0 
cííh,  MroKvne*  10 poner  f 

HRITEí*t300) 

FORMATE 1M* . ’Om  al  ancMc  an  aja  * s*  i»  tr*o«uién  , m 
RE AD I-, 400 > ANCHO* 

FQPrtATlFlO.D) 

IFÍ  ÚNCHOí  .LE.O.Oi'ñNCHOy-i  .o 
MLl  QOTOX Y i£t  I O , PONER  I 
HRITÉla.SQO) 

FfflWfttilH+p'duist*  Piltre  p**a  alte  í|>  o Piltra  pee*  bijDB  <£• 

ALAD  1*,4*>0}  TIPO 
*OC  FdRNRT  (F  10,03 

IF  ! ITIPOhNE. I .0) ,Í)R+íTIFO.NE*e.0TÍ  £30  TO  10 
CALL  HQTDItVte,  10,  PONER  j 
WR I TE  £ ■ f 700 1 

70O  FOfi*AmH+t‘ürcl*Á  d-t  filtre  ■ ’ h 

ftE  AD i ■ , BOa i ORDEN 
BO0  FQHfiñT  CF  10.01 

IF  (ÚRDEMELE, O. OI  QRDCN-l.Q 
00  EO  1-J¡,33 
M £0  3-1.33 

11- 36*1 

12- 34-1 

ie-34-j 

OJf-FLDAT £ 1-1  ) /CORTE* 

DT-FLQflTI J-l I/CORTEV 

&“CA,ea  < cmpl  m d « f dv  u 
IF  triLTito.Me.ai  theh 

IFIOAOEH.NÍ  „ | .O.AHIK.D.nE.0.0]  D-D**0ftt>EN 
K-0  ■ O 

IF  i F í LTRO . ED. 1 J K»|»0/f|.D  + 0,413*01 
*0l<F,LtBE0+5  íftMDj-  O *LT . ICO)  I-EíFE -0,347 

IF  t FfLTR0+EQ,4  i f-D 
ENDIF 

X -ORDEN 

JFlDat.  1,0)  60  TQ  22 

D*  í D— J . 0 1 -CORTE  S 

K-l -O  - ORDEN 

IF  í 0 , GT  .r  ANCHO*  ) EO  TQ  ££ 

Jt-1  ANCHQX-JJ)  /ANCHQk 
IFÍTIPÍJ.ED*  i,0|  I- 1 ,0  _ X 

IF  í ORDEN.  NE,t,0jjf.XaaORllCN 
aa  continué 

IF  CD-É0.0.0  .AND.  FILTRO. K.4)  X-TIPfl 


200 

300 

400 

10 

500 


ss 


LLENO  AHORA  LOS  FIXÉL'S 


IF  Í J I .01,03  1 SO  'TQ  40 
IF1  II . EQ . tó  ) 00  TQ  41 
AB  E 1 1 . j I í - ¡f 

41  S F 1 íe.Ea.-tjSi  GO  TO  4Ó 

ABC  IB, J1I*K 

40  IFÍ  J2.Eia.4i-  ) UO  TO  PO 
IFÍ  | I ,£0,43  i-  BO  TO  42 

apixi.jsi-x 

46  ¡FE  1?. £0,43»  90  TO  60 

AB(  12*321-* 

EO  CONTINUÉ 

REfllRN 
ENB 


su&RUjTiNE  transfon  1 1 nü  ice  i 

CÚHHON  /BÜFFÉft/  A 
COMPILES  AÍ*4,*41,AAí4Q0'0)  ,Nrtfn 
EOU1VALENCE  ÍA.AA1 
REAL  A¿A < *A , ¿4 i*RA i 409*1 ,99 1*4,441 
EOUIVALENtÉ  ( A , AA , RAA ,f¡A  I , C AAi  £049 1 , 99  1 
IF  ( iNOtCC.CO. L } THEN 
C TRANSFORMABA  DIRECTA  DE  FFT. 

PRINT  Triflirenidt  dirtcl»  da  f cunar  t FFT  |‘ 

SIGNO-* 1 .0 
K*1 

BO  100  1-1,44 

Ka3Ht 

DO  IOO  J*K,44,e 
AÍ  I,Jla-AC l,J) 

□O  IIO  1*1,44 

HALL  FA6TFFT (A( I , t 1 , SIQNU I 
CONTINUE 
DO  IEO  ¡f- I ,44 
DO  135  1*1,44 
WQRKl I 1-AEJ, I) 

CONTINUE 

CALI  FASTFFHHCMW. , SIGNO! 

DO  ia*  1-1,44 
M¡,  t l-W0Fit  í | ) 

cwit [mué 

COWT  | uC 

fc-t 

00  i 30  1-1,04 
K-3-K 

oo  130  j~+.,b4,a 
Al  I.JI—  Ai  1.31 
CONTINUE 

RETIKM 

ELSEIF  I ]NI}TCE.EO-,SE  THEN 
C TRANSFORMADA  INVERSA  DÉ  FtHJAlER. 

SIGNO-- I .O 
K*I 

DO  B00  1-1,44 
K-3-K 

DO  B00  J-S ,44 , B 

SCO  Ai  I , J )-"A( I , J ) 

DO  310  1-1 ,44 


IOO 

110 
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CALL  FASTFFT  t A 1 1 , I > I BNO t 

CONTINUE 

OO  £EO  J-í .¿4 

oo  agg 

rfDRl  ¿ I i »fl  (3,1  I 
CONTINUE 

CALÍ.  FASTFFT  ÍHOflF  .GIÜMO! 

BO  33A  I»U^ 

au»I  f-uqnM ! i 
continué 
CQNT INUE 

¡í-l 

DO  £30  1-1,44 

K-3-* 

DO  £30  J-|C,64,£ 

Al  t,  Jí— RU,J¡ 

CÜnT  |Hj£ 

PETU«N 

ENDIF 

END 


SU&NÜUTInC  FA9TFFT  I A,fll0NO¡ 
confie  i A í *4  >,¿1,  DOBLEN,  TE 

DATA  NH1  ,NV£ .N.h.PI/63,  35,  ¿***,3.  i 4 1 SígéÍÍ9SWW3¿ 

3- i 

OÚ  100  1-1  ,WI 

IF  r [ . GE .,  4 1 EO  TQ  101 

AUT-AÍ  I i 
Al IÍ-TE 
lOl  K-NVÉ 

líe  |F<r,,ÉK*4)  00  T(J  IC13 

kMc/e 
GO  Tq  IOS 
103  J-J+tC 
l<W  CONTINUE 

DO  UO  L-S  ,H 

L£-C*#i. 

LEI-LE^S 

u-li ,0,0*0) 

ANOUL  0-f*  I /FLÜAT  i LEI  ) 

ÜMIi|.EII"CNI'UfTOgg(«llWJ^OTl-&1GNO»ftr»<<  ANG<JLO>  > 

DO  1£<J  3—1  «LEI 
BQ  ¡13  1-4, N, LE 
ir-i+LEi 

te-ai iP>-o 

miFÍ-Aíll-TE 
Al I ¡-Al 1 J-TÉ 
113  CONTINUE 

l&>  CONTINUÉ 
1 10  CONTINUE 

DO  130  1*1, £.4 
ñül-A(  | )JQ.O 
130  CONTINUE 
END 


e i ■ ■ 
E£3 
5S* 

e®o 

@30 
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gUBftOtiTINE  IMOUTIOEC  ION.K.  É>UFfEft. CANAL  J 
CGhPl  E IC  BUFFERINl 
INTEOEfl  OFCiaN.CflNflJL 


Lftfi&Ü=44«N 
|F  lOPCION.EO.ll 

TT*N 

□a  IOD  1-1 ,L ASGO ,64 

lOO 

Kl*l 

RE  AD t CANAL  J i &UF  F ER I K I . K»fc 1 ,E£  J 

CONTINUE 

RETURN 

ELBEIF  ÍOFCION.EI0.£(  ¡MÉN 
L-Q  aoo  1-1  , LARDO, *4 
Kl-J 
T£-l  *4.3 

Uftt TE I CANAL J í BUTFSRTKl r K=KI.KE¡ 

COO 

CONTINUE 

FtETUFríJ 

CNOIF 

PRINT  ERROAí  hk»  •nlHt*  I»  o-pt  ion  qu* 

«■  Indice  Bfí  una  dm  la*  p«rwtro«: 

EW 
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NIVfUfl  MEN,  NftV  a O « 000 0 PQO E+OG  U3fi,D¡»0 

-32 
-31 
-Sil 
-2® 

-í! 

-a? 

-as 

-24 

-23 

-£2 

-21 

-20 

-13 

-ie 

-17 

-ie 

-i* 

-12 

-ja 


-10 

-t 

-e 


-a 

-i 

G 

1 

2 

3 

4 
3 
C 
7 

5 
3 


ia 

15 

14 


Fifi.  2.  Módulo  de  íu  tranxformedQ  de  FoUrier  de  la  Fig-  1. 
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HILELES  MIN,  HflX  1 tl.OOC!aD1íO£+flO  t + OQODQQ 
-32 
-31 
-30 
-29 
-20 
-27 
-2i 
-23 
-M 

-a 

-a? 

-ai 

-30 

-i? 

-ií 

-i? 

-Sé 

-13 

-H 

-15 

-11 

-10 

-5 

-a 


-* 

-5  4U¿U¿UJ, 


5 


£ 

7 

B 

9 

m 

11 

12 

13 

14 

15 
lí 
17 
IB 
1» 


ftg.  3r  Tcsf  rectangular. 
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Imagen,  Representación  espacial  de  un  objeto,  escena  u oirá  imagen. 
El  nivel  de  grií  I(x.y>  correspondiente  a la  coordenada  espacial  (K.y)  es,  en 
general,  proporcional  a la  energía  radíame. 

Nivel  de  grli.  Corresponde  al  número  o valor  asignado  a la  posición 
(x,y)  de  una  imagen.  Este  número  es  proporcional  a la  reflectarle  ia  o irans 
mitancia  de  una  pequeña  área  denominada  píxel. 

Intensidad  radiante.  De  un  punto  o un  objeto  es  una  medida  de  la 
potencia  radiante  por  ralcrcurradián  radiado  o reflejado  por  un  objeto.  La 
intensidad  radiante  es  una  función  de  la  naturaleza  del  objeto,  ángulo  de 
observación,  longitud  de  onda  espectral  y ancho  de  banda. 

tteflec  Lancia  o coeficiente  de  reflexión.  Es  el  cociente  entre  la  ener- 
gía por  unidad  de  tiempo  y por  unidad  de  área  (densidad  de  potencia  ra 
diante}  reflejada  por  un  objeto  y la  energía  por  unidad  de  tiempo  y por  uni- 
dad de  área  incidente  en  el  objeto. 

TransmJtancia  o coeficiente  de  irunsmil ancla.  Es  el  cociente  enlrc 
la  energía  por  unidad  de  tiempo  y unidad  de  área  (densidad  de  potencia 
radiante}  transmitida  por  un  objeto  y la  energía  por  unidad  de  tiempo  y uni- 
dad de  área  incidente  en  el  objeto, 

{.'mitra  s le.  El  contrasle  de  un  punto  frente  aun  fondo  puede  medirse 
por  1}  su  razón  de  contraste,  que  es  el  cociente  entre  el  valor  mayor  de 
la  transmiiancia  del  objeto  o del  fondo  y el  valor  menor  de  la  tnmsmitan- 
cia  del  objeto  o del  fundo,  y 2)  diferencia  de  Contraste,  que  es  la  diferencia 
entre  el  valor  mayor  de  la  transitó  lañe  i a del  objeto  o del  fondo  y el  valor 
menor  de  la  transmi  tañe  la  del  objeto  o del  fondo. 

Densidad.  La  densidad  en  un  punto  (x,y)  de  una  fotografía  es  una  me- 
dida de  la  capacidad  de  absorción  de  la  luz  de  La  plata  o del  tinte  deposi- 
tado en  cada  posición, 

Densil  órne  tro.  Dispositivo  utilizado  para  medir  el  valor  medio  de  la 
densidad  de  la  imagen  correspondiente  a una  pequeña  área  de  una  foto- 
grafía o de  una  transparencia.  Cuando  el  área  es  menor  que  unos  cente- 
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nares  de  mieras  cuadrada*,  el  dispositivo  recibe  el  nombre  de  micro- 
densitómetro. 

Resolución.  El  término  resolución  hace  referencia  a cómo  un  siste- 
ma. proceso,  componente  o material  puede  reproducir  un  objeto  aislado 
o discriminar  entre  objetos  próximos.  El  límite  de  resolución  o resolución 
espacial  viene  dada  por  la  dimensión  más  pequera  del  objeto  que  puede 
ser  discriminada. 

Poder  de  resolución.  El  poder  de  resolución  de  un  sistema  de  trata- 
miento  de  imágenes  es  una  medida  de  su  capacidad  de  discriminación  en- 
tre objetos  próximos.  El  poder  de  resolución  se  mide  en  pares  de  lineas 
por  milímetro,  es  decir,  el  máximo  número  de  lincas  y espacios  por  milí- 
metro que  pueden  reconocerse. 

Imagen  de  vídeo.  Es  una  imagen  que  presenta  un  formato  capaz  de 
ser  presentada  en  un  tubo  de  rayos  catódicos, 

Cuan  tífica  ció  n.  Proceso  mediante  el  cual  cada  nivel  de  gris  de  una 
imagen  fotograba,  etc.,  se  le  asigna  un  valor  entre  un  conjunto  finito  de 
niveles  de  gris. 

Flftni  espacial.  Transformación  de  una  imagen  utilizada  para  redu- 
cir el  ruido  o mejorar  algunas  características  de  la  imagen. 
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